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ABSTRACT 


Refinements  to  the  bow  and  stern  seal  subroutines  of  the 
XR-3  Loads  and  Motion  computer  program  are  made.   New  lift 
fan  performance  maps  are  developed  and  presented. "  Comparison 
studies  of  simulated  performance  of  the  original  XR-3  model 
with  the  old  and  new  fan  maps  are  made  under  varying 
sea  conditions.   Measured  roll  and  pitch  transient 
behavior  of  the  XR-3  test  craft  in  calm  water  is  compared 
with  the  simulated  response  of  both  the  original  and  the  new 
model.   In  addition,  simulated  performances  of  the  original 
and  new  models  in  irregular  ahead  seas  are  compared. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  speed  limitation  of  the  conventional  displacement 
vessel  has  been  recognized  for  many  years.   In  an  effort  to 
effect  a  quantum  Increase  In  surface  vessel  speed  the  United 
States  Navy  has  been  developing  various  craft  whose  principal 
means  of  support  Is  other  than  hydrostatic  lift. 

One  type  of  craft  currently  receiving  much  attention  Is 
the  surface  effect  ship  or  craft.   In  this  category  are  the 
air  cushion  vehicles  and  captured  air  bubble  craft  (CAB). 
The  general  nature  of  these  craft  and  their  construction  is 
well  presented  by  Robert  L.  Trillo  in  Ref.  1.   In  these 
craft  either  all  or  a  major  portion  of  its  support  is  derived 
from  a  pressure  differential  between  the  atmosphere  and  a 
plenum  chamber  which  is  open  at  the  bottom.   The  speed 
advantages  gained  by  this  type  of  vessel  are  derived  from 
two  principal  characteristics.   It  does  not  waste  energy 
displacing  a  large  volume  of  v;ater  and  the  frictional  forces 
betvjeen  the  hull  and  water  are  greatly  reduced  by  the  craft 
riding  on  a  cushion  of  air  with  only  a  minimum  of  the  craft 
structure  actually  in  contact  with  the  water. 

Vessels  of  this  type  are  generally  separated  into  two 
categories  -  "Air  Cushion  Vehicles"  whose  weight  is  entirely 
supported  by  the  pressure  differential  In  the  plenum  chamber 
and  "Captured  Air  Bubble  Craft"  whose  weight  is  partially 
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supported  by  sldewalls  which  extend  Into  the  water.   It  Is 
the  second  type  of  vessel  In  which  the  United  States  Navy 
has  become  primarily  interested  for  ocean  going  applications. 

The  term  Captured  Air  Bubble  (CAB)  is  possibly  a  misnomer 
since  in  this  type  of  craft  the  air  does  leak  out  and  has  to 
be  continually  replenished  by  supply  fans.   The  leakage 
though  is  relatively  slow  when  it  is  compared  to  the  air 
cushion  vehicle  which  has  continuous  air  gap  around  the 
periphery  of  the  craft  and  hence  the  fans  that  replenish  the 
air  in  the  plenum  chamber  have  to  be  larger  and  more  powerful 
than  those  of  a  CAB  of  similar  size. 

This  thesis  is  concerned  with  one  aspect  of  the  develop- 
ment of  the  Captured  Air  Bubble  vehicle  and  that  is  the 
sim.ulation  of  such  craft  utilizing  a  digital  computer, 
specifically  the  Loads  and  Motion  Program  developed  by 
OceanicSj  Inc. 

The  basic  rigid  body  analysis  and  spatial  relationships 
of  the  Loads  and  Motion  program  are  well  documented  in 
Ref.  2.   Also  covered  are  the  principal  dynamic  and  static 
approximations  used  in  developing  the  equations  of  motion 
for  the  craft  in  its  six  degrees  of  freedom. 

A  copy  of  this  program  was  delivered  to  the  Naval 
Postgraduate  School  in  October  of  1972  which  had  been 
specifically  tailored  for  the  Bell  (lOOB)  100  ton  craft. 
Leo  and  Boncal  in  Ref.  3  succeeded  in  converting  this  program 
to  represent  the  XR-3  test  craft. 
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B.   OBJECTIVE 

The  purpose  of  this  thesis  was  to  take  a  second  more 
detailed  look  Into  all  of  the  aspects  of  seal  modeling  and 
the  fan  performance  map  with  the  Idea  of  possibly  making 
improvements  In  these  representations  before  testing  the 
computer  model  of  the  XR-3  under  dynamic  conditions. 
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II.   FAN  PERFORMANCE  MAP  ANALYSIS  AND  COMPARISON 

A.   INTRODUCTION 

As  the  surface  effect  ship  is  predominately  dependent 
on  the  plenum  lifting  force  for  its  performance,  the 
characteristics  of  the  lifting  fans  are  important  if  valid 
studies  of  the  craft  in  question  are  to  be  made. 

Leo  and  Boncal  reported  in  Ref.  3  that  the  steady  state 
operating  conditions  of  the  XR-3  were  obtained  by  increasing 
the  leakage  area  of  the  craft  until  the  proper  equilibrium 
pressure  was  maintained  in  the  main  plenum.   This  area  was 
twenty-six  square  feet.   As  the  cross-sectional  area  of  the 
plenum  is  less  than  twenty  square  feet  this  leakage  area  was 
out  of  proportion  to  properly  represent  the  leakage  area  of 
the  XR-3. 

While  operating  above  hump  speed  the  XR-3  test  craft  has 
been  observed  to  have  a  gap  underneath  the  sternseal  of 
from  four  to  five  inches.   Slight  leakage  does  occur  at 
the  bowseal  as  reported  by  Ref.  3  but  the  total  leakage  area 
is  estimated  to  be  not  over  four  square  feet. 

To  appreciate  the  effect  this  leakage  area  has  on  computed 
craft  performance  an  understanding  is  required  of  the  method 
by  v;hlch  the  pressure  is  calculated  in  the  plenum. 

When  the  program  is  executed  an  estimation  of  the  initial 
operating  conditions  must  be  supplied  by  the  user.   Besides 
velocity  and  thrust,  Initial  conditions  for  plenum  pressure, 
displacement  and  pitch  angle  must  be  supplied.   From  the 
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displacement,  pitch  angle  and  velocity  for  calm  water 
conditions,  a  plenum  volume  is  computed.   Using  the  leakage 
areas  supplied,  exit  flow  rates  are  computed  for  the  bow 
and  stern  seals.   The  fan  subroutine  is  then  called  using 
the  pressure  in  the  plenum,  bow  and  stern  seals  as  primary 
arguments.   Using  these  pressures  supply  flow  rates  are 
established  for. the  fans  that  supply  air  to  the  plenum, 
and  both  seals.   As  the  seals  discharge  into  the  main  plenum 
the  sum  of  these  fan  flow  rates  establish  the  total  influx 
rate.   The  net  flow  being  the  difference  between  the  influx 
and  efflux  rates  establishes  the  new  plenum  pressure  for 
the  next  calculation  interval  based  upon  the  adiabatic 
perfect  gas  law.   This  new  pressure  creates  a  new  plenum 
lifting  force  which  is  then  balanced  by  the  defining  differ- 
ential equations  of  motion  resulting  in  a  new  operating 
condition  of  craft  attitude  and  displacement. 

In  studying  craft  motion  in  calm  water  the  leakage  area 
is  of  little  or  no  consequence  since  no  additional  leakage 
areas  are  normally  developed  beneath  the  sidewalls  or  seals, 
and  the  rate  of  change  of  plenum  volume  is  relatively  slow. 
On  the  other  hand,  when  the  program  is  executed  in  simulated 
sea  states  where  gap  areas  are  generated  beneath  the  seals 
and  sidewalls  and  plenum  volume  does  change  quite  rapidly 
the  leakage  areas  and  fan  maps  become  of  prime  importance. 

It  must  be  emphasized  though  that  it  is  not  the  leakage 
area  Itself  that  affects  craft  performance.   It  is  rather 
the  relative  size  of  the  fixed  leakage  area  with  respect  to 


the  variable  leakage  areas  generated  around  the  periphery 
of  the  craft  by  wave  motion  and  craft  attitude  and  position 
which  Is  Important. 

B.   FAN  MAP  ANALYSIS  AND  DERIVATION 

Isolation  of  the  factors  which  caused  the  large 
discrepancy  between  observed  steady  state  gap  area  and 
calculated  area  required  an  examination  of  both  the  efflux 
and  Influx  equations  for  possible  errors. 

The  equation  used  to  establish  leakage  flow  rates  under 
the  seals  and  the  sldewalls  are  the  same.   This  equation 
is  given  below 


Q  =  C^A^(2Ap/p)^/^ 


It  is  the  general  flow  equation  for  a  sharp  edged  orifice 
in  turbulent  flow  conditions.   As  the  orifice  coefficient, 
C„,  is  the  only  variable  in  the  equation  that  could  be 
dependent  on  craft  geometry  and  its  range  was  limited  between 
unity  and  six  tenths;  the  efflux  equations  were  eliminated 
as  the  source  of  the  error.   The  only  area  left  for  possible 
error  was  the  fan  maps  themselves. 

The  fan  maps  in  question  are  arrays  which  give  fan  flow 
rates  in  cubic  feet  per  second  as  a  function  of  plenum  or 
seal  pressure.   They  are  entered  as  data  in  "INCON"  as 
characteristics  of  the  specific  craft  being  modeled. 
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It  was  found  that  the  fan  maps  constructed  by  Leo  and 
Boncal  to  represent  the  fans  installed  In  the  XR-3  were  based 
upon  Incorrect  information.   Actual  fan  performance  data  for 
the  fans  currently  installed  is  presented^ in  Figure  1,  and 
physical  characteristics  in  Figure  2. 

The  fan  performance  map  given  by  the  manufacturer,  Joy 
Manufacturing  Company,  is  based  upon  flow  that  is  not  impeded 
by  air  ducts.   It  is  therefore  optimistic  for  any  system 
that  contains  long  ducts  as  these  ducts  will  modify 
considerably  the  fan  performance  by  creating  ducting  losses. 

The  best  way  to  obtain  the  actual  fan  performance  maps 
in  a  duct  system  would  be  to  actually  measure  the  flow  rates 
as  the  pressure  in  the  receiving  reservoir  (in  this  case  the 
plenum)  is  varied.   As  this  method  was  impractical  in  the 
case  of  the  XR-3  another  method  of  estimating  the  modified 
fan  performance  had  to  be  devised. 

The  method  outlined  below  was  utilized  and  although  it 
probably  does  not  give  exact  results  the  fan  maps  thus 
achieved  should  adequately  represent  the  actual  fan  performance 
of  the  XR-3. 


The  fan  data  was  provided  by  Associate  Professor  of 
Aerology  D.M.  Layton  v/ho  is  in  charge  of  the  actual  XR-3 
test  program  at  the  Naval  Postgraduate  School.   He  has 
verified  that  the  fans  Installed  are  Joy  Model  AVR87-55D2058 
for  which  the  data  given  in  Figures  1  and  2  are  correct. 
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FIGURE   1 


Performance  and  Fan  Map  Curves  for  the  Joy  Ax i vane  Fan 
(Model  AVR  87-55D2058) 
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The  first  step  used  In  arriving  at  the  new  fan  maps  was 
to  establish  a  realistic  leakage  area.   Once  this  leakage 
area  was  established  a  steady  state  efflux  rate  was 
established.   Because  in  steady  state  operating  conditions 
the  influx  and  efflux  rates  must  be  equal  the  flow  rate 
supplied  by  the  lift  fans  in  total  was  established.   It  was 
at  this  point  that  a  second  assumption  was  made.   This 
assumption  was  that  there  were  greater  losses  in  the  ducting 
supplying  air  to  the  bow  seal  and  main  plenum  chamber  than 
there  were  in  the  ducting  supplying  the  rear  seal.   This 
assumption  v;as  based  upon  the  fact  that  the  ducting  to  the 
bowseal  and  main  plenum  had  to  carry  the  air  forward 
approximately  sixteen  feet  to  their  discharge  orifices 
whereas  the  stern  seal  ducting  was  only  a  few  feet  long  with 
only  one  curve  in  it  to  Impede  flow. 

As  a  result  the  fraction  of  total  flow  assigned  to  the 
stern  seal  was  greater  on  an  individual  fan  basis  than  was 
the  air  supplied  by  the  bov;  seal  or  main  plenum  fans. 
Another  assumption  was  that  the  same  map  could  adequately 
represent  both  the  bow  seal  supply  fans  and  those  supplying 
the  main  chamber.   This  last  assumption  was  based  on  similar 
duct  lengths  and  construction. 

Now  knowing  the  individual  fan  flow  rates  and  the  actual 
plenum  pressure  a  pressure  drop  in  the  duct  for  steady  state 
conditions  was  established.   This  pressure  drop  was  the 
difference  between  the  known  plenum  operating  pressure 
(2i].8i<  lbs/ft^)  and  the  indicated  fan  discharge  pressure 
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for  that  flow  rate.   This  pressure  drop  was  then  scaled 
for  different  flow  rates  based  upon  the  fact  that  pressure 
losses  vary  with  the  square  of  the  flow  rates.   In  this 
manner  new  fan  maps  were  computed  for  the  bow  seal  and 
main  plenum,  and  the  stern  sea.   These  maps  are  presented 
in  Figures  3  and  ^. 

The  last  step  was  to  go  back  to  the  leakage  area  and 
make  the  slight  adjustments  needed  to  give  the  required 

plenum  steady  state  pressure  of  2^1.8^  lbs/ft  .   The  final 

2 
leakage  areas  established  were  .08  ft   for  the  bow  seal 

2 
and  3.79  ft  . for  the  stern  seal. 

In  the  original  fan  subroutine  the  horsepower  generated 
by  the  fans  is  computed  based  upon  the  influx  rate  and  the 
pressure  differential  between  the  atmosphere  and  the  plenum. 
As  the  actual  work  done  per  unit  time  by  the  fans  this 
figure  can  be  misinterpreted  as  the  actual  power  needed 
by  the  fans  to  generate  printed  work  rate.   This  would  only 
hold  true  if  the  fans  operated  with  one  hundred  percent 
efficiency  and  there  were  no  ducting  losses.   In  actuality 
fans  of  this  nature  v;hich  supply  relatively  large  volumes 
of  air  at  lovj  pressure  rarely  operate  at  over  seventy  percent 
efficiency  and  can  achieve  that  only  over  an  extremely 
narrow  operating  range.   These  physical  characteristics  of 
fans  \-iere   well  presented  by  R.D.  Moyer  in  19'^6  in  Ref.  ^. 

For  the  fans  that  are  used  in  the  XR-3  the  Joy  Manufac- 
turing Company  has  provided  the  fan's  horsepower  requirements 
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FIGURE  3 

Bow  Seal  and  Plenum  Chamber  Fan  Maps^  adjusted  for  ducting  losses 
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FIGURE  H 

Stern  Seal  Fan  Map,  adjusted  to  account  for  ducting  loss 
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on  the  same  performance  graph  as  the  fan  maps.   This  power 
requirement  Is  plotted  for  the  entire  operating  range  of 
the  fans  and  Is  fairly  consistent  at  about  2  horsepower. 
This  horsepower  requirement  does  not  change  under  the  new 
fan  maps  since  the  fan  map  discharge  pressures  at  the  fan 
exhaust  remains  the  same.   This  horsepower  map  was  entered 
In  the  fan  subroutine  as  a  table  of  look  up  and  actual  power 
requirements  were  computed  from  it  using  flow  rate  as  the 
entering  argument. 

In  the  COLFIL  subroutine  (the  routine  that  prints  out 
all  data  and  draws  all  graphs)  the  fan  power  label  was 
changed  to  "fan  power  (ideal)"  to  avoid  possible  confusion. 
In  the  fan  subroutine  the  print  option  was  modified  to  give 
the  actual  power  required  as  well  as  the  ideal  power.  Also 
the  fan  efficiency  is  computed  and  is  given  as  an  output  of 
this  subroutine. 

C.   SIMULATED  PERFORMANCE  COMPARISONS  BETWEEN 
ORIGINAL  AND  NEW  FAN  MAPS 

As  stated  earlier  it  was  believed  that  changing  the 
leakage  areas  and  the  fan  maps  would  have  a  significant 
effect  only  when  the  computer  model  was  run  under  conditions 
where  the  plenum  volume  underwent  comparatively  rapid 
fluctuations  and/or  the  leakage  areas  around  the  periphery 
of  the  hull  varied. 

These  assumptions  v;ere  confirmed  by  simulated  runs  under 
varying  sea  conditions.   In  all  these  comparison  tests  the 
two  XR-3  models  used  were  identical  except  that  the  original 
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fan  maps  were  replaced  by  the  newly  developed  ones  In  one 
model.   Both  simulated  craft  were  run  using  the  original 
seals  as  developed  In  Ref.  3.   Craft  weight  was  5813  lbs, 
center  of  gravity  9.57  feet  from  the  stern,  speed  was  held 
constant  at  2  0  knots.   All  runs  were  made  into  ahead  seas. 
The  sea  conditions  chosen  were  irregular  sea  state  1, 
irregular  sea  state  2,  regular  seas  1  and  regular  seas  2. 
The  latter  two  were  so  named  because  the  amplitudes  and 
frequencies  chosen  for  these  regular  seas  runs  were  the  major 
component  of  irregular  sea  state  1  and  irregular  sea  state  2 
respectively.   The  actual  amplitudes  and  frequencies  for 
these  runs  are  tabulated  on  page  27. 

In  the  simulation  runs  into  regular  seas  1  dramatic 
differences  were  noted  in  the  fluctuation  of  plenum  pressure 
and  heave  accelerations.   With  the  original  fan  maps, 
variations  in  plenum  pressure  were  limited  to  about  2.2  lbs/ft' 
whereas  with  the  new  maps  the  variations  were  over  6.9  lbs/ft 
This  represents  an  increase  in  fluctuations  of  over  315 
percent.   Because  the  plenum  normally  provides  over  seventy 
percent  of  the  supporting  force  for  the  craft  similarly 
large  variations  were  noted  in  heave  acceleration  both  at 
the  bov/  and  at  the  center  of  gravity.   Pitch  motion,  though, 
did  not  change  much  being  less  than  .1  degree  in  both  cases. 

Results  for  the  runs  in  other  sea  conditions  are 
tabulated  on  page  23.   The  differences  noted  for  the  runs 
in  irregular  sea  state  2  and  regular  seas  2   v.'ere  not  as  great. 
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This  has  been  attributed  to  the  craft's  ability  to  conform 
its  attitude  to  the  wave  slope  in  the  case  of  the  longer 
wave  lengths  associated  with  greater  sea  states.   This  ability 
to  conform  then  reduces  the  high  frequency  volumetric  changes 
in  the  plenum  chamber,  thereby  reducing  pressure  fluctuations. 
Figures  5  through  12  are  the  results  of  running  both  models 
in  the  simulated  sea  state  one.   By  visually  comparing  the 
corresponding  variances  of  the  two  models  it  is  readily 
apparent  that  the  fan  maps  do  significantly  change  the 
performance  characteristics  of  the  model.   By  examining  the 
table  on  page  28  it  is  also  apparent  that  one  must  be  careful 
in  predicting  behavior  as  the  variation  in  measured  variables 
does  not  vary  directly  with  the  sea  conditions.   The  new  fan 
maps  have  in  effect  defined  a  similar  but  different  craft. 
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TABLE  I.   SEA  STATE  WAVE  COMPONENTS 


Run  1   Regular  Seas  1   Heading  l80 

Wave  Components      Amplitude  (ft) 
1  .0598 


Frequency  (rad/sec) 
4.1526 


Run  2    Regular  Seas  2   Heading  iBO 

Wave  Components      Amplitude  (ft) 
1  .2721 


Frequency  (rad/sec) 
1.8447 


Run  3   Irregular  Sea  State  1   Heading  l80 

Wave  Components      Amplitude  (ft)    Frequency  (rad/sec) 


1 

.0029 

2.553 

2 

.0324 

3.256 

3 

.0599 

4.1526 

4 

.0558 

5.2961 

5 

.0401 

6.7545 

6 

.0262 

8.6145 

7 

.0165 

10.9867 

8 

.0102 

14.0122 

Rim  4   Irregular  Sea  State  2   Heading  l80 

Wave  Components      Amplitude  (ft)    Frequency  (rad/sec) 


1 

.0168 

1.1786 

2 

.1464 

1.4745 

3 

.2721 

1.8447 

4 

.2668 

2.3077 

5 

.2052 

2.8871 

6 

.1410 

3.6118 

7 

.0928 

4.5185 

8 

.0600 

5.6528 
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TABLE  II   SIMULATED  XR-3  PERFORMANCE  WITH  DIFFERENT 

FAN  MAPS 


Comparison  of  XR-3  Loads  and  Motion  Computer  Program 
performance  with  original  for  maps  with  XR-3  with  new 
fan  maps  are  given  below: 


RUN 

QUANTITY 

XR-3 
(original) 

XR-3 
(new) 

1. 

Regular 
Seas  1 

Plenum  Pressure 

Pitch  Angle 

Draft 

C.G.  Acceleration 

24.8±1.09  PSF 
3±.05  deg. 
6.25±.05  in. 
±.04  g 

24.8±3.45  PSF 
2.7±.05  deg. 
5.9±.l  ^. 
±.15  g 

2. 

Ix'regular 
Sea  State  1 

Plenum  Pressure 

Pitch  Angle 

Draft 

C.G.  Acceleration 

24.8±2.0  PSF 
3.2±.l  deg. 
6.5±.05  in. 
±.07  g 

24.8±15  PSF 
2.5±.2  deg. 
5.75±.l  in. 
±.4  g 

3. 

Regular 
Seas  2 

Plenum  Pressure 

Pitch  Angle 

Draft 

C.G.  Acceleration 

23.75±5.75  PSF 
5.1±4.75  deg. 
7.25±2.3  in. 
±.23  g 

24.8±7  PSF 
4.2±3.8  deg. 
6.0±2.9  in: 
±.27  g 

4. 

Irregular 
Sea  State  2 

Plenum  Pi -es  sure 

Pitch  Angle 

Draft 

C.G.   Acceleration 

5  to  33  PSF 
.8±2.2  deg. 
8±6  in. 
±.46  g 

-5  to  40  PSF 
.8±2.5  deg. 
10±7  in. 
±:7  g 
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FIGURE    5 

Computer  Simulation  Plenum  Pressure  verses  Time.  Speed  20  knots 
Original  Fan  Maps,  Sea  State  1     .     . 
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FIGURE    6 

Computer  Simulation  Pitch  Angle  verses  Time,  Speed  20  knots 
Original  Fan  Maps,  Sea  State  1  . 
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FIGURE  7 

Computer  Simulation  C.G.  Draft  verses  Time,  Speed  20  knots 
Original  Fan  Maps,  Sea  State  1 
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FIGURE  8 

Computer  Simulation  C.G.  acceleration  verses  Time.  Speed  20  knots 
Original  Fan  Maps>  Sea  State  1 
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FIGURE  9 

Computer  Simulation  Plenum  Pressure  verses  Time^  Speed  20  knots 
New  Fan  Maps,  Sea  State  1 
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FIGURE    10 

Computer  Simulation  Pitch  Angle  verses  Time>  Speed  20  knots 
New  Fan  Maps^  Sea  State  1 
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FIGURE    11 

Computer  Simulation  C.G.  Draft  verses  Timc.,  Speed  20  knots 
New  Fan  Haps,  Sea  State  1 
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FIGURE  12   ■  ■ 

Computer  Simulation  C.G.  Acceleration  verses  Time.  Speed  20  knots 
New  Fan  Maps,  Sea  State  1 
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III.   SEAL  MODEL  REFINEMENTS 

A.   INTRODUCTION  AND  BASIC  MODELING  TECHNIQUES 

Through  tests  conducted  aboard  the  XR-3  by  the  NFS  test 
group  headed  by  Associated  Professor  D.M.  Layton  of  the 
Aeronautics  Department  It  has  been  demonstrated  that  the 
forces  generated  by  the  seals  have  an  appreciable  effect  on 
the  overall  performance  of  the  test  craft.   By  varying  the 
lower  travel  limit  of  the  bow  seal  and  contou'ring  its  shape 
the  steady  state  pitch  angle  can  be  changed  as  well  as  the 
thrust  versus  velocity  characteristics.   Although  not  as 
pronounced,  changing  the  support  cable  length  of  the  stern 
seal  also  affected  craft  trim  conditions. 

Because  of  the  demonstrated  steady  stats  pitch  sensitivity 
of  the  XR-3  to  seal  position  limits  as  set  by  the  support 
cables  it  was  decided  that  more  detailed  modeling 'of  the 
seals  v/as  needed  before  the  overall  response  of  the  L&M 
computer  model  would  reflect  the  behavior  of  the  XR-3  under 
varying  conditions  of  support  cable  length. 

The  original  XR-3  seal  modeling  as  presented  in  Ref.  3 
was  patterned  after  the  rear  seal  in  the  L&K  computer  model 
of  the  Bell  (lOOB)  100  ton  CAB.   This  seal  is  of  the  "soft" 
pneumatic  type  and  its  representation  is  well  covered  by 
Kaplan,  Bentson  and  Sargent  in  Ref.  2. 

The  pneumatic  "soft"  seal  model  is  based  on  the  assumption 
that  this  type  of  seal  is  unstlffened  and  ncn-rlgid.   The 
mass  of  the  seal  was  considered  so  insignificantly  small  as 
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compared  to  the  forces  generated  by  the  seal  when  It  is  In 
contact  with  the  water  as  to  have  negligible  effect  on  overall 
craft  dynamics.   The  lifting  forces  generated  by  the  seals 
were  calculated  by  multiplying  the  wetted  contact  area  of 
the  seals  by  the  pressure  differential  between  the  seal  and 
the  atmosphere  for  the  bow  seal  and  the  plenum  chamber  and 
seal  in  the  case  of  the  stern  seal. 

Drag  forces  on  the  seals  were  considered  to  be  due  only 
to  skin  friction  between  the  water  and  the  seal.   The  equation 
used  to  compute  these  forces  was 


F      =  -  -^  C^.  A    , 
^seal      2   di   seal 


where  p  is  the  density  of  water,  u  the  velocity  of  the  craft 

and  A    -,  is  the  wetted  surface  of  the  seal.   C,„  is  the 
seal  df 

drag  coefficient  and  is  calculated  using  the  Schultz-Griinow 
formula 


^df  "  -^^^^   (logio^^  "  -^O?)^'^^ 


where  Re  is  the  Reynolds  number  computed  based  on  the  average 
wetted  length  of  the  seal  and  the  Reynolds  number  is 

calculated  by : 

T3   _  u  X  (wetted  length) 
ne  —  F 

1.28x10"^ 


where  1.28x10"^  ft  /sec  is  the  kinematic  viscosity  of  the 
water. 
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The  actual  construction  of  the  seals  Is  well  described 
In  Ref.  3.   Each  seal  consists  of  two  pneumatic  lobes 
separated  by  a  common  membrane.   The  membrane  has  limber 
holes  in  it  to  equalize  the  pressure  between  the  lobes.   The 
lower  surface  of  the  bottom  lobe  is  shaped  by  twelve  steel 
battens  that  act  as  stiffeners.   These  battens  are  spaced 
uniformly  across  the  bottom  surface  at  right  angles  to  the 
hinge  line.   It  is  this  surface  that  comes  in  contact  with 
the  water,  and  will  hereafter  be  referred  to  as  the  "leading 
face"  of  the  seal. 

The  forward  edge  of  the  leading  faces  is  hinged  near  the 
top  of  the  plenum  chamber  allowing  the  seal  to  pivot  about 
this  line  when  waves  are  encountered. 

The  pressure  in  the  seals  is  maintained  by  axial  fans. 
The  air  is  delivered  by  ducts  that  discharge  into  the  top 
of  the  upper  lobe.   The  seals  in  turn  vent  to  the  main 
plenum  by  means  of  ducting.   Back  pressure  in  the  vent 
ducting  can  be  controlled  by  butterfly  valves  that  can  be 
manually  adjusted.   Normally  the  seals  are  operated  so  that 
they  maintain  about  1-2  PSF  more  pressure  than  the  plenum. 

Downward  travel  limits  of  the  seals  are  controlled  by 
adjustable  support  cables  attached  to  the  rear  edge  of  each 
steel  batten  at  the  trailing  edge  of  the  seal's  leading  face. 
In  addition  a  second  support  cable  is  attached  to  each 
stiffener  approximately  midway  between  the  hinge  line  and 
the  trailing  edge.   Figure  13  shows  a  cutaway  of  the  seals, 
only  two  cables  are  dravm  in  for  clarity. 
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Both  the  bow  and  the  stern  seal  extend  the  width  of 
the  plenum  chamber  and  are  Identical  in  construction.   The 
leading  faces  are  approximately  ^9  inches  long  from  hinge 
line  to  trailing  edge.   The  perpendicular  distance  from 
sldewall  keels  to  the  hinge  line  is  22.5  inches. 

As  stated  before  the  original  modeling  of  the  seals  did 
not  take  into  account  any  of  the  seal's  ability  to  flex  or 
pivot  on  its  hinge.   In  addition  no  attempt  was  made  to 
model  the  support  cables  and  their  ability  to  affect  craft 
attitude.   In  the  new  modeling  the  support  cables  have  been 
incorporated.   In  the  stern  seal  modeling  only  the  rear 
support  cable  is  modeled  since  this  seal  rarely  comes  in 
contact  with  the  water.   In  the  bow  seal  subroutine  both  the 
rear  support  cables  and  the  center  support  cables  have  been 
included  since  this  seal  does  contact  the  water.   Before  any 
attempt  at  revising  the  seal  subroutines  could  be  undertaken 
it  was  necessary  first  to  study  the  action  of  the  seals  while 
the  craft  was  actually  underway. 

The  semi-rigid  construction  of  the  seal  dictated  a  highly 
complex  distributed  parameter  model  of  the  seals  be  developed 
if  the  seals  were  to  be  fully  modeled  without  any  approxima- 
tions.  In  addition  to  the  fact  that  practically  all  of  the 
physical  characteristic  parameters  of  the  seal  such  as  spring 
constants  for  the  steel  stiffeners  and  fabric  flexibility 
v;ere  unknovm  such  a  model  would  require  an  inordinate  amount 
of  calculation  time  in  the  computer.   Since  extra  computation 
time  was  to  be  avoided  simplifying  assumptions  had  to  be  made. 
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To  be  accurate  these  assumptions  or  approximations  had  to 
be  based  on  actual  observations  of  the  seal  or  on  previously 
made  approximations  whose  effect  has  already  been  Investigated. 

In  Ref.  2  the  approximation  of  eliminating  seal  dynamics 
and  the  setting  of  seal  pressure  at  a  constant  level  above 
plenum  pressure  has  been  thoroughly  Investigated  and  was 
Incorporated  In  the  original  model  of  soft  pneumatic  seals. 
These  approximations  have  been  continued  In  the  present  seal 
modeling. 

The  other  modeling  approximations  made  were  based  upon 
actual  observations  made  of  the  seals  while  the  XR-3  was  in 
motion.   These  observations  are  described  below: 

1.  For  small  v;ater  perturbations  (approximately  3  inches) 
the  bcv;  seal  flexed  and  appeared  to  pivot  slightly  about  the 
center  support  cable  causing  the  angle  at  the  hinge  line 

to  slightly  increase,  but  almost  insignificantly. 

2.  For  longer  perturbations  the  seal  both  flexed  as 
described  above  and  then  began  to  pivot  on  the  hinge.   After 
the  seal  rose  about  six  Inches  no  additional  flexing  took 
place. 

3.  The  bov;  seal  maintained  contact  with  the  water  at  all 
times  and  never  gaped  because  of  an  inability  to  follow  wave 
motion. 

^.      The  bow  seals  penetration  of  the  water  was  practically 
independent  of  water  motion  in  relation  to  the  craft,  but 
appeared  to  be  inversely  related  to  the  craft's  velocity. 
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When  the  XR-3  was  at  rest  penetration  was  about  five  to  six 
Inches  and  when  traveling  at  20  knots  appeared  to  be 
approximately  1  Inch. 

5.   The  stern  seal  rode  out  of  the  water  while  the  craft 
operated  above  "hump"  speed.   The  height  of  the  air  gap  thus 
created  appeared  to  remain  almost  constant  at  four  to  five 
inches.   Water  rarely  contacted  the  seal. 

Using  the  above  described  observations  and  the  preceding 
simplifying  approximations  concerning  seal  dynamics  and 
pressures  both  the  bow  and  stern  seal  were  modified  to 
reflect  the  action  of  the  actual  seals. 

The  approximation  given  in  Ref.  2  of  equating  the  wetted 
length  of  the  seal  with  the  water  height  above  the  sidewall 
keel  at  the  seal's  trailing  edge  location  was  eliminated. 
In  its  stead  a  two-dimensional  "table  of  look  up"  was 
developed  which  gave  wetted  length  as  a  function  of  the 
immersion  depth  of  the  seal's  trailing  edge  in  the  water  and 
the  seal's  trailing  edge  deflection  height  above  the  keel  at 
the  trailing  edge  location.   In  other  words  if  the  seal  was 
deflected  upwards  because  of  wave  motion  or  pitching  motion 
of  the  craft  the  deflection  height  of  the  trailing  edge 
would  increase  and  the  seal's  leading  face  would  become  more 
parallel  to  the  water's  surface  thereby  Increasing  the  wetted 
area  of  the  seal. 
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The  "table  of  lookup"  was  constructed  using  the 
information  obtained  from  a  graphical  deflection  model 
similar  to  the  one  presented  in  Figure  14.   In  this  model 
the  leading  face  of  the  seal  Is  represented  in  varying 
amounts  of  deflection  height  from  full  down  (point  A)  to 
full  up  against  the  seal  framework  (point  B).   Full  down 
is  zero  deflection  and  is  where  the  trailing  edge  of  the 
seal  is  even  with  the  bottom  of  the  sldewall  keel,  full  up 
is  maximum  deflection.   This  deflection  height  was  divided 
into  segments  of  one  inch  for  ease  of  computation.   Points 
were  marked  upon  each  leading  face  line  which  represented 
the  surface  of  the  water  on  the  leading  face  under  varying 
depths  of  trailing  edge  immersion.   The  distance  along  the 
leading  face  to  the  trailing  edge  from  these  points  then 
represented  the  length  of  the  leading  face  that  was  immersed 
in  the  v/ater  or  the  "wetted  length".   The  points  for  similar 
immersion  depths  v;ere  connected  for  the  varying  seal  deflec- 
tion heights  and  these  construction  lines  are  labeled  with 
their  corresponding  immersion  depths  in  inches  along  the 
zero  deflection  height  leading  face.   The  construction  line 
corresponding  to  an  immersion  depth  of  zero  Inches  is,  of 
course,  congruent  v/lth  the  trailing  edge  of  the  seal  since 
if  the  seal  did  not  penetrate  the  water  to  at  least  some 
degree  there  could  be  no  "wetted  length"  unless  the  leading 
face  was  exactly  parallel  to  the  water's  surface,  which 
never  occurs.   These  construction  lines  v;ere  drawn  for 
Immersion  depths  of  the  trailing  edge  of  0  to  5  Inches. 
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The  actual  Immersion  depth  of  the  trailing  edge  is 
calculated  based  upon  the  observation  that  at  rest  the 
immersion  depth  was  approximately  5.5  inches  and  while 
traveling  at  20  knots  was  estimated  at  a  little  less  than 
one  inch.  Since  the  hydrodynamic  force  on  the  seal  is 
proportional  to  the  square  of  the  craft's  velocity  the 
immersion  depth  was  assumed  to  be  approximately  inversely 
proportional  to  some  factor  of  the  square  of  the  speed. 

Using  5.5  inches  as  the  known  immersion  depth  at  zero 
zpeed  and  .75  inches  at  20  knots,  the  following  empirical 
formula  was  derived  which  gives  reasonable  results 


D(ln.)  =  5.5/(1  +  ^)2 


where  D(ln. )  is  the  depth  in  inches  and  u  is  the  velocity 
in  feet  per  second. 

Using  the  calculated  immersion  depth  and  the ' seals 
calculated  deflection  height  the  table  is  entered  to  give 
the  wetted  length.   This  table  of  look  up  is  presented 
graphically  in  Figure  15.   The  non-linear  nature  of  the 
curves  is  readily  apparent.   An  Important  aspect  of  this 
table  apparent  in  Figure  15  is  that  within  about  the  first 
ten  Inches  of  seal  deflection  the  v;etted  length  of  the  seal 
remains  fairly  constant,  especially  if  the  craft  is  operating 
at  speeds  where  the  immersion  depth  is  approximately  one  inch. 
This  figure  also  shows  that  the  maximum  wetted  length  is 
equal  to  the  length  of  the  leading  face  (just  over  48  Inches). 
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As  the  seal  deflects  and  becomes  more  horizontal  the  wetted 
surface  area  becomes  greater  given  a  constant  Immersion 
depth. 

This  same  "table  of  look  up"  is  provided  in  the  stern 
seal  modeling  but  because  the  seal  does  not  ride  in  the  water 
an  assumed  "depth"  had  to  be  used  as  an  entering  argument. 

The  assumption  was  made  that  since  the  stern  seal  does 
not  ride  on  the  water  its  lifting  force  has  to  be  the  result 
of  some  form  of  air  pressure.   This  force  could  only  be 
generated  by  the  impingement  of  air  particles  on  the  seal's 
leading  face  as  they  rushed  to  escape  underneath  the  seal. 

At  this  point  an  assumption  was  made  that  the  effective 
wetted  area  occurred  where  the  kinetic  energy  in  the  air  mass 
was  approximately  one  quarter  that  of  the  air  escaping  under- 
neath the  seal.   As  the  air  in  the  plenum  is  at  such  a 
pressure  that  it  can  be  considered  incompressible  this  would 
occur  where  the  cross-sectional  area  was  approximately  two 
times  the  cross-sectional  area  beneath  the  seal.   As  the 
width  of  the  plenum  chamber  is  constant  and  the  gap  height 
is  approximately  4.5  inches  at  all  times  the  point  on  the 
seal  where  the  cross-sectional  area  is  two  times  the  gap 
area  is  the  point  along  the  leading  face  where  the  height 
from  the  seal  to  the  water  is  two  times  the  gap  height  or 
about  9  Inches.   Subtracting  the  gap  height  gives  an 
effective  "phantom"  immersion  depth  of  4.5  inches;  this  was 
rounded  to  five  inches. 
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Lift  forces  generated  by  the  stern  seal  were  considered 
to  be  the  product  of  the  seal  area  found  by  using  the  above 
phantom  immersion  depth  and  the  seal  deflection  height,  and 
the  pressure  differential  between  the  stern  seal  and  plenum 
chamber. 

The  above  lifting  force  calculation  has  only  minimal 
effect  on  simulated  performance  of  the  XR-3  since  the 
pressure  between  the  seal  and  the  plenum  chamber  is  normally 
only  one  to  two  PSF.   The  drag  forces  on  the  stern  seal  are 
calculated  in  the  same  manner  as  they  are  in  the  bow  seal 
subroutine  but  are  subsequently  set  to  zero  since  the  seal 
rarely  contacts  the  water. 

B,   BOW  SEAL  MODELING 

Contained  herein  is  a  detailed  analysis  of  the  actual 
bow  seal  subroutine.   It  is  coded  in  standard  FORTRAN.   By 
reviewing  the  actual  subroutine  the  actual  approximations 
used  in  the  routine  may  be  understood.   A  full  listing  of 
the  BOV/SL  subroutine  is  contained  in  Appendix  A  along  with 
the  required  entries  in  the  data  input  statements. 

(1)  DO  1  J=l,ll 

(2)  GAP(J)  =  0.0 

(3)  ELSKI(J)  =  0.0 

(4)  WETLEN(J)  =0.0 

(5)  ELSKID(J)  =  0.0 

(6)  1  CONTINUE 

The  above  statements  (1  through  6)  zero  all  variable  arrays 
that  will  be  used  every  time  BOWSL  is  entered.   GAF(J)  is 
the  height  of  the  air  gap  at  each  station  along  the  bow  seal, 
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ELSKI(J)  Is  the  water  height  above  the  keel  at  the  trailing 
edge  of  the  seal.   WETLEN(J)  Is  the  wetted  surface  length 
at  each  station.   ELSKID(J)  Is  the  deflection  height  of  the 
seal  trailing  edge  measured  from  the  botton  of  the  keel  of 
the  craft. 

(7)  ALBS  =0.0 

(8)  FX  =  0.0 

(9)  FY  =  0.0 

(10)  FZ  =  0.0 

(11)  FK  =  0.0 

(12)  FM  =  0.0 

(13)  FN  =  0.0 

The  preceding  seven  statements  Initialize  the  single  valued 
variables  that  are  computed  each  time  BOWSL  is  entered. 
ALBS  is  leakage  area  of  the  bow  seal.   FX  is  the  forces 
parallel  to  the  craft's  axis  of  symmetry  (centerline) . 
FY  is  the  horizontal  forces  perpendicular  to  the  FX  forces. 
FZ  is  the  vertical  forces.   FK,  FM,  FN  are  the  moments  about 
the  Z  axis,  Y  axis  and  X  axis  respectively. 

(I'D  DELPBG  =  PBS  -  PB 

(15)  IF  ( DELPBG. LT. 0.0 )DELPBG=0.0 

Equation  l4  computes  the  pressure  differential  between 

the  plenum  chamber  (PB)  and  the  seal  (PBS).   This  differential 

is  actually  fixed  in  INCON  as  an  entering  argument.   Equation 

15  assures  that  the  seal  differential  pressure  is  never 

negative.   This  is  reasonable  since  the  seals  vent  to  the 

chamber. 
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(16)  ARGO  =  ELMAXB/CORLEN 

(17)  ANGO  =  ARSIN(ARGO) 

(18)  XI  =  XBS+ZBS*THETA-CORLEN*COS(ANGO) 

The  three  equations  above  (I6-I8)  establish  the  moment 
arm  in  the  X  direction  from  the  craft  center  of  gravity 
to  the  trailing  edge  of  the  bow  seal.   ELMAXB  Is  the  rear 
cable  length.   CORLEN  Is  the  chord  length  of  the  bow  seal. 
XBS  is  the  distance  from  the  center  of  gravity  to  the  bow. 
ZBS  Is  the  vertical  distance  of  the  bow  hinge  line  below 
the  center  of  gravity.   Theta  is  the  pitch  angle  of  the 
craft. 

(19)  Zl  =  -Z-ABS+XBS«THETA-ELMAXB«COS(THETA) 

Equation  (19)  establishes  the  moment  arm  in  the  -Z 
direction  to  the  water's  surface  through  which  the  drag 
forces  act  to  create  a  moment  about  the  Y  axis. 

(20)  DPTHFT  =  (5. 5/(1. +(U/25.))**2. 0/12.0 

Equation  (20)  calculates  the  immersion  depth  of  the 
seal  based  upon  craft  velocity,  u. 

(21)  IP(CENCAB.GT.l. 1875 )CENCAB=1. 1875 

Logical  statement  (21)  sets  the  maximum  length  of  the 
center  support  cable  to  3  inches  more  than  the  maximum  rear 
cable  length,  ELMAXB. 
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(22)  DO  3  K=1,N 

3  CONTINUE 

The  above  loop  calculates  for  each  seal  partition  point 
the  water  height  above  the  keel  ELSKl(K),  the  wetted  length, 
WETLEN(K),  and  the  seal  deflection  height  ELSKID(K)  in 
preparation  to  calculating  the  forces  and  moments. 

(23)  DPFT(K)  =  DPTHFT 

Equation  (23)  initializes  the  immersion  depth  of  the 
seal. 

(24)  ELSK1(K)=(ETA(3,K)-DTABX(K)«(XX(3,K)-X1)-Z1) 

+yY ( 3 ,K) *PHI+XLXPWV*WATSLP 

Equation  (2^1)  establishes  the  water  height  above  (or 
below)  the  keel  at  the  bow.   ETA(3,K)  is  the  height  of  the 
water  at  the  center  of  gravity  above  mean.   DTABX(K)  is 
the  correction  factor  for  moving  this  height  away  from  the 
C.G.  (Center  of  Gravity).   (XX(3,K)-X1)  is  the  distance  the 
correction  factor  has  to  act  on.   Zl  is  the  height  of  the 
keel  in  relation  to  the  C.G.   YY(3,K)  is  the  displacement 
of  the  partition  from  the  craft  center  line.   PHI  is  the 
roll  angle.   XLXPNV*WATSLP  is  the  correction  factor  for  the 
water  slope  caused  by  the  plenum  pressure  acting  for  a 
finite  period. 
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(25)     IF(ELSKI(K).GT.HINGHT)  ELSKI(K)=HINGHT 

Equation  (25)  limits  the  water  height  to  the  hinge 
height  of  the  bow  seal. 


(26)     IF((HINGHT-ELSKI(K)+DPFT(K)).GE.ELMAXB) 
DPFT(K)=ELMAXB-HINGHT+ELSKI(K) 


Equation  (26)  adjusts  the  immersion  depth  of  the  seal 
so  that  the  seal  can  go  no  lower  than  that  permitted  by  the 
rear  cable  ELMAXB. 

(27)  IF(DPFT.LT.O.O)  DPFT=0.0 

DPFT  is  used  in  computing  DPIN  which  is  an  entering 
argument  into  the  wet  length  table  BWSL.   Therefore 
equation  (27)  limits  this  argument  to  a  minimum  value  of 
zero. 

(28)  ELSKID(K)=(ELSKI(K)-DPFT(K) )*12 .0 

Equation  (28)  calculates  the  deflection  height  of  the 
seal  (ELSKID(K))  in  Inches  based  upon  water  height,  ELSKI(K) 
and  the  immersion  depth  DPFT(K). 

(29)  IF( (HINGHT-ELSKID(K)/12. 0) .GE. ELMAXB) 

ELSKID(K) = (HINGHT-ELMAXB) «12 . 0 

Equation  (29)  limits  the  downward  travel  of  the  bow  seal 
to  ELMAXB  the  rear  cable  length  and  adjusts  the  deflection 
height  ELSKID(K)  to  reflect  this  lower  limit. 
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(30)  ■  DPIN=DPFT(K)*12.0 

(31)  MM=DPIN 

(32)  MM1=MM+1 

(33)  MM2=MM1+1 
(3'»)     DINC=DPIN-MM 

The  equations  above,  (30)  through  (3^),  establish 
parameters  for  entering  the  bow  seal  wetted  length  array, 
BOWSL.  DPIN  Is  the  immersion  depth  In  Inches,  MMl  Is  the 
lower  limit  table  calling  Integer,  MM2  Is  the  upper  table 
calling  Integer  and  DINC  Is  the  incremental  difference. 

(35)  GAP(K)=-ELSKI(K)+(HINGHT-ELMAXB) 

Equation  (35)  computes  the  air  gap  height  based  on 
water  height  ELSKI(K)  which  is  positive  upwards  and  the 
limit  of  the  seal  travel  (HINGHT-ELMAXB) . 

(36)  IP(GAP.LT.O.O)  GAP=0.0 

Equation  (36)  limits  the  gap  to  positive  values. 

(37)  IF(ELSKID(K) .GE.0.0)  GO  TO  2 

(38)  WETLEN(K)=ELSKI(K) 

(39)  GO  TO  3 

Equations  (37)  through  (39)  set  the  wet  length  equal  to 
the  water  height  If  it  is  negative.   This  computation  is 
used  in  case  one  seal  partition  point  has  a  negative  value 
while  the  others  on  either  side  are  positive. 
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(40)  MM3=ELSKID(K) 

(111)  MMi^=Mr43+l 

(k2)  MM5=MM4+1 

(i»3)  DLING=ELSKID(K)-MM3 

Equations  (kO)    through  (^13)  establish  the  final  two 
entering  arguments  for  retrieving  the  four  points  needed 
from  BOWSL  to  compute  the  wetted  length.   ELSKID(K)  Is  the 
deflection  height  of  the  seal.   MM4  Is  the  lower  entering 
Integer  based  upon  the  deflection  height,  MM5  Is  the  upper 
integer.   DLINC  Is  the  Incremental  difference  used  for  linear 
interpolation. 

.     ■  (4i|)  BWSL1=BWSL(MM1,MM'4) 

(45)  BWSL2=BWSL(MM1,MM5) 

(46)  BWSL3=BWSL(MM2,MM4) 

(47)  BWSL4=BV/SL(MM2,MM5) 

(48)  BWSLA1=(BWSL2-BWSL1)*DLINC+BWSL1 

(49)  BV/SLA2=(BVJSL4-BVJSL3)*DLINC+BWSL3 

(50)  WETLEN(K)=( (BWSLA2-BWSLA1) *DINC+BWSLA1)/12 . 0 

Equations  (44)  through  (50)  identify  the  four  points  of 
the  BWSL  array  which  bracket  the  actual  value  of  wetted 
length.   Linear  interpolation  is  used  twice  to  find  the 
average  based  on  deflection  height  and  once  to  find  the 
final  wetted  length.   The  array  BWSL  was  calculated  in  inches; 
thus,  to  convert  to  feet  the  final  average  is  divided  by 
12.0. 

(51)  DO  8  J=1,N 

8  CONTINUE 
This  loop  makes  all  the  force  calculations  based  upon 
wetted  length  and  seal  position  calculated  previously.   In 


55 


addition  It  calculates  the  leakage  area  based  upon  the  steady 
state  leakage  area  and  the  variable  gap  heights. 

(52)  WETLAV=(WETLEN(J+l)+WETLEN(J))/2.0 

(53)  IF(WETLAV.LE. .001)00  TO  6 
(5^)     DPFTAV=(DPFT(J+l)+DPFT(J))/2.0 

(55)  ELSKIA=(ELSKI(J+l)+ELSKI(J))/2.0 

(56)  ELSKDA=(ELSKID(J+1)+ELSKID(J) )/2.0 

Equations  (52)  and  (5^)  through  (56)  average  the 
wetted  length,  immersion  depth,  water  height  and  deflection 
height  respectively,  for  each  seal  segment  based  on  the 
values  calculated  for  the  stations  on  either  side. 
Equation  (53)  routes  the  program  around  the  remaining  force 
calculations  if  the  wetted  length  of  the  seal  is  zero  or 
insignificant  for  calculation. 

(57)  SEALHT=niNGHT-ELSKDA 

Equation  (57)  is  the  calculation  of  the  distance  between 
the  top  of  the  seal's  supporting  frame  and  the  bottom  of 
the  trailing  edge  of  the  seal. 

(58)  .       DIFF=2.0  5fCENCAB-(SEALHT+0.5) 

(59)  IF(DIFF.GT.0.5)    DIFF=0.5 

The  DIFF  equation  (58)  calculates  the  spring  distance 
of  the  steel  battens  based  upon  the  center  support  cable 
length  CENCAB  and  the  position  of  the  seal,  SEALHT. 
Equation  (59)  limits  the  spring  movement  to  .5  feet  or 
six  inches,  the  point  where  rigid  body  analysis  comes  into 
full  effect. 
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(60)  ARMlB(J)=Xl+WETLAV/2.0 

(61)  ARM2B ( J )  =  ZS-ELSKI A+DPFTAV/2 . 0 

Equations  (60)  and  (6l)  calculate  the  moment  arms  in 
the  X  and  Z  direction  respectively  based  upon  the  length  to 
the  trailing  edge  of  the  bow  seal,  XI,  and  half  the  wetted 
length,  WETLAV/2.0,  for  ARMIB.   ARM2B  Is  based  upon  the 
height  of  the  center  of  gravity  ZS ,  the  water  height, 
ELSKIA  and  half  the  Immersion  depth  of  the  seal,  DPFTAV/2.0. 

(62)  IF(DIFF.GE.0.25)  GO  TO  i| 

(63)  DFBS(J)=-DELP-DELYBS*WETLAV 

(64)  GO  TO  5 

Equations  (62)  through  (64)  calculate  the  vertical  force 
acting  on  the  segment  If  the  deflection  difference  Is  less 
than  .25  feet  or  three  Inches.   If  the  deflection  difference 
Is  greater  than  .25  or  3  Inches  the  vertical  force  equation 
has  a  spring  constant  term  In  It  acting  upon  the  moment 
summations  about  the  hinge  line. 

(65)  4  FORLEN=XBF-WETLAV 

(66)  IF(FORLEN.EQ.OO)  GO  TO  5 

(67)  ARGV/=  CHINGHT-ELSKIA  )  /FORLEN 

(68)  IF(ARGV/.GT.1.0)ARGW=1.0 

(69)  ANGV/-ARSIN(ARGV;) 

(70)  FORCOS=COS(ANGW) 

(71)  GO  TO  6 

(72)  5  FORCOS=0.0 

Equations  (65)  through  (72)  calculate  the  length  of  the 
seal  not  in  contact  with  the  water,  FORLEN,  and  the  angle 
which  this  length  makes  with  the  water  ANGW.   From  the  angle 
the  cosine  of  ANGV/  is  calculated,  FORCOS,  for  future  use 
in  calculating  the  vertical  moments  about  the  hinge  line. 
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(73)         DFBS(J)=-DELP*DELYBS*WETLAV-DELP*FORLEN*DELYBS 

*FORCOS  * ( FORLEN  *0 . 5  *FORCOS ) / 
FORLEN*FORCOS+WETLAV/2.0) )* 
((DIFF-0.25)*4.0) 


The  above  equation  calculates  the  vertical  force  acting 
upon  the  bow  seal  by  summing  vertical  moments  about  the 
hinge  line  and  equating  It  to  the  moment  created  by  one 
vertical  force  vector  acting  on  the  midpoint  of  the  wetted 
•surface.   The  term  "  (DIFF-0.  25)  *4.  0"  Introduces  a  spring 
distance  of  three  inches  through  which  the  vertical  force 
transitions  from  being  only  based  upon  the  product  of  the 
wetted  area  and  the  pressure  differential  to  one  based  upon 
equalizing  moments  about  the  hinge  line.   After  the  seal 
bends  this  three  Inches  It  Is  assumed  to  be  completely  rigid 
and  therefore  If  the  seal  Is  not  undergoing  accelerations, 
the  moments  must  balance  about  the  hinge  line. 

(7^)  ARG=0.5*RHO*U*U«WETLAV«DELYBS 

(75)  RESKI=U*WETLAV/ENU 

(76)  CDTSKI  =  0.il2  7/(ALOG10(RESKI)-0.il0  7)-'^*2.6  4 

(77)  TSKIB(J)=-ARG«CDTSKI 

The  above  four  equations  calculate  the  drag  forces  of 
the  seal  based  upon  the  previously  discussed  Schult z-Grunow 
formula.   RESKI  is  the  Reynolds  number ,CDTSKI  Is  the  drag 
coefficient  and  TSKIB(J)  is  the  drag  force  on  the  segment. 
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7  CONTINUE 

(78)  FX=FX+TSK1B(J) 

(79)  FZ=FZ+DFBS(J) 

(80)  FK=FK+DFBS(J)*YAVGB(J) 

( 81 )  FM=FM-DPBS ( J ) *ARM1B ( J ) +TSKB ( J ) « ARM2B ( J ) 

(82)  FN=FN-TSK1B(J)*YAVGB(J) 

(83)  ALBS=ALBS+GAP(J)+GAP(J+l))*DELYBS/2,0 

Equations  (78)  through  (83)  calculate  the  forces  and 
moments  for  each  segment  and  sums  them  to  give  total  forces 
and  moments.   YAVGB(J)  is  the  average  distance  off  centerline 
of  the  individual  segment. 

(8H)  8  CONTINUE 

(85)  ALBS=ALBS+BLEAK 

(86)  SQFAC=SQRT(2.*ABS(PBAR)/RH0INF) 

(87)  QL=CFBS*ALBS«SQFAC*SIGN(1. ,PBAR) 

The  last  equations  calculate  the  flow  rate  of  air  based 
upon  the  leakage  area  ALBS  and  the  pressure  differential 
between  the  plenum  chamber  and  the  atmosphere.   BLEAK  is 
the  seal  leakage  entered  in  INCON  to  account  for  slight 
leakage  where  the  seals  do  not  quite  maintain  contact  with 
the  sidewall. 

The  remaining  statements  of  subroutine  BOWSL  give  the 
following  values  as  printed  output  if  the  bov;  seal  print 
switch  is  on,  GAP,  WETLEN,  FX,  FY,  FZ ,  FK,  FM,  and  FN. 

C.   STERN  SEAL  MODELING 

The  stern  seal  subroutine  is  covered  in  the  same  manner 
as  the  bov;  seal  subroutine.   The  complete  listing  of  the 
stern  seal  subroutine  is  contained  in  Appendix  B  along  with 
instructions  on  the  proper  entries  to  be  used  on  the  data 
cards. 
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(1)  DO  1  J=l,ll 

(2)  GAP(J)=0.0 

(3)  ELSKI(J)=0.0 
(1)      ELSKIL(J)=0.0 

(5)  AIRLEN(J)=0.0 

(6)  1  CONTINUE 

The  above  loop  and  equations  zero  all  variable  arrays 
that  will  change  each  time  subroutine  STNSL  is  entered  as 
was  the  case  in  the  bow  seal  subroutine. 

(7)  EFDEP  =6.0 

Equation  (7)  establishes  the  effective  depth  of  the  seal 
at  five  inches.  This  figure  is  used  as  an  entering  argument 
in  the  stern  seal  contour  mapping  array  CTNSL. 


(8) 

ALSS=0.0 

(9) 

FX=0.0 

(10) 

FZ=0.0 

(11) 

FK=0.0 

(12) 

FM=0.0 

(13) 

FN=0.0 

(14) 

AGAP1=0.0 

(15) 

AGAP2=0.0 

(16) 

AGAPA1=0.0 

The  nine  equations  above  zero  all  single  valued  variables 
that  are  computed  each  time  STNSL  is  entered. 

(17)  DELP=PSS-PB 

(18)  IF(DELP.LT.O.O)  DELP=0.0 

The  two  statements  above  fix  the  pressure  differential 
between  the  plenum  (PB)  and  the  stern  seal  (PSS).   Equation 
(17)  is  first  used  to  arrive  at  the  difference  between 
these  two  quantities  and  the  logical  IF  statement  is  used 
to  insure  all  values  are  non-negative.   As  in  the  bow  seal 
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subroutine  this  is  logical  since  the  air  supplied  to  the 
stem  seal  vents  into  the  plenum  chamber. 

(19)  PBAR=PB-PINF 

This  equation  establishes  the  pressure  differential 
between  the  plenum  (PB)  and  the  atmosphere.   It  is  used  in 
computing  leakage  flow  rates. 

(20)  SINDIF=SINTH-COSTH*THETA 

(21)  COSDIF=COSTH+SINTH*THETA 

The  above  equations  (20)  and  (21)  establish  the  angle 
of  the  XR-3  with  respect  to  the  horizontal  and  the  vertical, 
respectively. 

(22)  X1=XSS+ZSS*THETA-XLF*SINDIF 

(23)  Z1=(-Z-ZSS+XSS*THETA-ELMAXS«C0S(THETA)) 

The  above  two  equations  establish  lengths  in  the  X  and  Z 
directions  which  are  used  to  calculate  the  water  height  above 
the  keel  at  the  stern  seal  location. 

(24)  N=NSTA(4) 

(25)  DO  2  K=1,N 
2  CONTINUE 

The  above  loop  is  used  to  establish  the  gap  at  each 
station  and  the  phantom  "wetted  length"  of  the  seal. 
NSTA(4)  is  the  number  of  stations  chosen  for  the  seal  and 
can  be  up  to  eleven,  thereby  defining  ten  seal  segments. 
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(26)  ELSKI(K)=(ETA(4,K)-DETADX(K)»(XX(4,K)-X1) 

-Z1)+YY(4,K)*PHI 

(27)  IF(ELSKI(K) .GT.HINGHT)    ELSKI(K)=HINGHT 

Equation  (26)  computes  the  height  of  the  water  above 
the  keel  at  the  stern  seal  location.   ETA(^,K)  Is  a  length 
correction  factor  that  is  dependent  on  seal  shape.   Since 
in  the  case  of  the  XR-3  the  seal  is  flat  across  the  hinge 
line  this  term  is  zero.   DETADX(K)  is  a  wave  correction 
factor  that  is  multiplied  by  the  seal  distance  from  the 
center  of  gravity  (XX(4,K)-X1)  to  give  the  wave  height  at 
the  stern  seal  location.   The  logical  IF  statement  limits 
the  water  height  to  the  top  of  the  seal. 

(28)  ELSKIL(K)=ELSKI(K)+GPS 

(29)  IF(ELSKIL(K) .GT.HINGHT)  ELSKIL(K)=HINGHT 

(30)  IF(ELSKIL(K) .LT. (HINGHT-ELMAXS) ) 

ELSKIL(K)=HINGHT-ELMAXS 

ELSKIL(K)  is  the  height  of  the  trailing  edge  of  the 
stern  seal  above  the  keel.   It  is  measured  in  feet.   GPS  is 
the  air  gap  height  that  was  observed  to  be  constant  under 
normal  running  conditions.   It  is  .4644  feet.   The  logical 
IF  statement  (29)  establishes  the  upper  limit  of, seal 
deflection  where  it  is  stopped  by  the  seal  supporting 
framev/ork.   The  second  logical  IF  statement  establishes  the 
lower  limit  of  seal  deflection  height. 

(31)  GAP(K)  =  -ELSKI(K)  +  (HINGHT-ELr.^AXS) 

(32)  IF(GAP(K).LT.O.O)GAP(K)-0.0 

The  equation  (31)  and  logical  IF  statement  (32)  establish 
the  gap  height  below  the  lower  stop  limit  set  by  the  support 
cable  length  ELMAXS. 
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(33)  MM1=ELSKIL(K)*12.0 

(3^)  MM2-MM1 

(35)  MM3=MM2+1 

(36)  DLINC=ELSKIL(K)*12.0-MM1 

(37)  STNSL1=CTNSL(MM,MM2) 

(38)  STNSL2  =  CTNSL(MM,Mr43) 

(39)  AIRLEN(K)=((STNSL2-STNSL1)*DLINC+STNSL1)/12.0 

Equations  (33)  through  (35)  establish  the  entering 
arguments  of  the  seal  contour  array  CTNSL  In  the  same  manner 
as  was  done  in  the  bow  seal  subroutine.   AIRLEN(K)  Is  the 
phantom  "wetted  length"  and  Is  given  the  name  AIRLEN 
because  the  seal  is  supported  by  air. 

(40)  N=NSTA(J4)-1 

(41)  DO  5  J=1,N 

5  CONTINUE  . 

The  above  loop  accomplishes  the  actual  force  and 
moments  of  each  seal  segment  and  totals  the  forces  and 
m.oments.   It  also  computes  the  average  gap  height  based  upon 
the  gap  heights  of  the  individual  stations. 

( 42 )  ELSKI A= (ELSKI ( J+1 ) +ELSKI ( J ) ) /2  .  0 

(43)  ELSKLA=(ELSKIL(J+l)+ELSKlL(J))/2.0 

(44)  AIRLAV=(AIRLEN(J+l)+AIREN(J))/2.0 

The  three  above  equations  average  the  water  height  above 
the  keel,  ELSKIAjthe  average  seal  deflection  height  ELSKLA, 
and  the  average  phantom  wetted  length  AIRLAV. 

(45)  AGAP=ELSKLA-ELSKIA 

(46)  AGAP1=AGAP 

Equation  (45)  establishes  the  average  height  of  the  air 
gap  at  the  seal  segment.   AGAPl  is  a  variable  used  in  further 
calculations. 
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(^7)     IP(AGAP.LT.GPS)AGAP=GPS 
(i\8)  IF(AGAP1.GT.GPS)AGAP1=GPS 

The  above  logical  IF  statements  establish  the  lower 
limit  for  AGAP  at  .HSHh   feet  at  the  maximum  value  of  AGAPl 
at  .4644  feet. 

(49)  ARMlS=XX(4,J)+ELSKIA/2 

(50)  ARM2S=ZS-ELSK1A 

The  above  equations  (49)  and  (50)  calculate  the  moment 
arms  through  which  the  vertical  lifting  force  acts  to  create 
the  major  component  of  the  FM  moment  about  the  Y  axis;  the 
second  equation  creates  the  moment  arm  with  respect  to  the 
C.G.  through  which  the  drag  forces  act. 

(51)  DFSS(J)=-DELP*DELYSS«AIRLAV/(GPS/AGAP)**2.0 

Equation  (51)  calculates  the  vertical  force  acting  upon 
the  seal  and  with  the  term  (GPS/AGAP) **2 . 0  allows  for  this 
force  to  decrease  with  the  square  of  the  ratio  of  the  actual 
gap  v;ldth  AGAP  and  the  nominal  width  GPS.   This  is  based  on 
the  assumption  that  the  air  exiting  beneath  the  stern  seal 
supports  the  stern  seal  and  that  the  lift  forces  generated 
by  the  exiting  air  are  Inversely  proportional  to  their 
velocity. 
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(52)  ARG=.5RH0*U*U*AIRLAV*DELYSS 

(53)  RESKI=U*AIRLAV/ENU 

(5^)  CDTSKI=.42  7/(ALOG10(RESKI)-.if07)**2.6i| 

(55)  TSKIS(J)=-ARG*CDTSKI 

(56)  TSKIS(J)=0.0 

Equations  (52)  through  (55)  calculate  the  theoretical 
drag  forces  based  upon  the  average  phantom  wetted  length 
AIRLAV  as  if  the  seal  was  actually  in  contact  with  the  water. 
The  Schultz  Grunow  formula  is  used  as  in  the  bow  seal  routine 
to  calculate  the  drag  forces  based  upon  the  Reynold's  number 
which  is  based  on  AIRLAV  the  average  "wetted"  length. 
Equation  (56)  removes  this  drag  calculation  since  the  seal 
is  presumed  to  ride  on  a  cushion  of  air  and  is  hence 
negligible. 

(57)  FX=PX+TSKIS(J) 

(58)  FZ  =  FZ-fDFSS(J) 

(59)  FK-FK+DFSS(J)»YAVGS(J) 

( 60 )  PM=FM-DFSS ( J ) *ARM1S ( J ) +TSKIS( J ) »ARM2S ( J ) 

(61)  FN=FN-TSKIS(J)«YAVGS(J) 

The  five  equations  above  calculate  the  forces  and 
moments  of  each  segm.ent  and  sums  them  for  all  the  segments 
to  give  the  total  seal  forces  and  mioments. 

(62)  ALSS=ALSS+(GAP(J)+GAP(J+l))*DELYSS/2.0 

(63)  AGAP2=AGAP2+AGAP1 

(64)  AGAPA1=AGAP2/J 

The  above  equation  (62)  calculates  the  total  gap  area 
ALSS  that  is  due  to  the  seal  reaching  the  downward  limit 
of  its  travel  and  the  gap  thus  increasing  past  its  nominal 
value  of  .4644  feet.   Equations  (63)  and  (64)  calculate  the 
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average  gap  height  including  the  nominal  .^6^^  height  and 
allows  the  nominal  gap  to  close  if  the  seal  is  slammed 
against  its  upper  limit,  the  upper  seal  support  framework. 

(65)  ALSS=ALSS+ALEAK*(AGAPA1/GPS) 

Equation  (65)  calculates  the  leakage  area  under  the  seal 
ALSS  is  the  leakage  area  due  to  the  gap  width  increasing 
beyond  the  nominal  .H6Hi\   feet.   The  term  ALEAK  is  the 
nominal  leakage  area  of  3.79  square  feet  normally  observed. 
Multiplying  ALEAK  by  AGAPAl/GPS  allows  this  nominal  leakage 
area  to  decrease  if  the  seal  is  limited  in  its  upward  travel, 

(66)  SQPAC=SQRT(2.»ABS(PBAR/RH0INF) 

(67)  QL=CFSS*ALSS*SQFAC»SIGN(1. ,PBAR) 

Equations  (66)  and  (67)  calculate  the  flow  rate  based 
upon  the  leakage  area  and  the  pressure  differential  between 
the  plenum  and  the  atmosphere.  The  equation  assumes  orifice 
type  leakage. 

The  remaining  statements  give  print  outputs  of  the  terms 
GAP,  AIRLEN,  FX,  FY,  FZ ,  FK ,  FM  and  FN  if  the  stern  seal 
print  sv;itch  is  on. 
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IV.   DATA  ANALYSIS  OF  ACTUAL  CRAFT  TRANSIENT 
MOTION  AND  COMPUTER  SIMULATION 

A.   INTRODUCTION 

In  order  to  verify  the  modeling  changes  made  In  the 
seals,  data  was  collected  on  the  XR-3  test  craft  for 
conditions  of  induced  transient  motion  in  pitch  and  roll. 
To  induce  roll  action  the  test  craft  was  sallied  by 
shifting  body  weight  laterally  near  the  C.G.  until  maximum 
oscillations  built  up,  then  the  forcing  action  was  stopped. 
For  pitch  oscillations  the  chase  boat  was  veered  across 
the  path  of  the  XR-3  to  produce  a  wave  front  with  its  wake. 
By  having  the  XR-3  pass  over  the  wake  it  was  intended  to 
induce  a  transitory  pitch  oscillation. 

These  attempts  to  produce  a  transient  response  cannot 
be  viev;ed  as  entirely  successful.   The  principal  reason  for 
the  limited  success  was  that  it  was  extremely  difficult  to 
induce  a  large  enough  perturbation  in  either  pitch  or  roll 
so  that  the  normal  background  oscillations  could  be  ignored. 

Figures  l6,  17  and  l8  are  graphs  of  the  roll  and  pitch 
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motions   of  the  XR-3  when  making  its  twenty  knot  calibration. 

The  magnitude  and  frequency  of  the  steady  state  pitch 


Figure  17  depicts  the  actual  pitch  angle  of  the  test 
craft.   Figure  l8  is  the  same  data  with  the  steady  state 
average  pitch  angle  removed.   This  was  done  on  all  data  to 
facilitate  the  measurement  of  the  pitch  oscillations  super- 
imposed on  the  steady  state  value  of  pitch  angle. 
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FIGURE  16 

Roll  Motion  verses  Time^  Velocity  20.^29  knots 
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FIGURE    17 

Pitch  verses  Time,  Velocity  20.^29  knots 
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FIGURE  18 

Pitch  Hotion  verses  Time,  Velocity  20.^29  knots 
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oscillation  shown  In  Figure  l8  is  approximately  .08-  degrees 
and  .295  Hz  respectively.   The  steady  state  magnitude  of 
the  roll  motion  shown  in  Figure  16  was  approximately  .075 
degrees  and  its  frequency  was  approximately  the  same  frequency 
of  the  pitch  motion,  .295  Hz. 

These  steady  state  oscillations  appear  quite  small  until 
It  is  realized  that  the  maximum  roll  perturbation  that  was 
able  to  be  induced  during  any  run  was  only  about  .^   degrees. 
Lower  ratios,  between  induced  and  background  oscillations 
were  the  norm  for  both  roll  and  pitch  induced  motion.   The 
background  oscillations  in  pitch  and  roll  therefore  could 
not  be  Ignored  when   analyzing  data,  and  deleterlously 
affected  the  measurement  of  decay  time  constants. 

As  stated  before  the  attempt  was  made  to  introduce 
separately  pitch  and  roll  motion  into  the  test  craft  and 
then  record  the  transient  response  after  the  forcing  action 
had  passed  or  was  stopped.   After  studying  the  processed 
data  It  was  evident  that  it  had  been  nearly  impossible  to 
affect  transients  in  one  variable  (either  pitch  or  roll) 
without  introducing  at  least  some  oscillations  in  the  other 
variable.   In  introducing  pitch  oscillations  roll  motion  was 
greatly  affected.   This  was  caused  because  the  wave  front 
caused  by  the  chase  boat's  wake  always  struck  the  test  craft 
at  an  oblique  angle  and  not  perpendicular  to  the  axis  of 
symmetry  of  the  craft.   Only  during  those  runs  where  the 
XR-3  was  sallied  was  some  modicum  of  success  achieved  In 
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Isolating  roll  transient  motion.   The  one  run  Included  In 
this  report  where  sallying  was  used  to  Induce  roll  transients 
Is  the  one  conducted  at  11.79  knots.   During  the  other  tv/o 
runs  Included  (15.93  and  19.^5  knots)  the  wake  crossing 
method  was  used  to  force  oscillations.   This  resulted  In 
runs  In  which  both  roll  and  pitch  transients  could  be 
studied  at  the  same  time. 

In  making  comparison  between  the  recorded  data  and  that 
of  the  computer  simulations  both  the  original  XR-3  model  as 
presented  In  Ref.  3  and  the  new  model  which  incorporates 
the  new  fan  maps,  new  bow  and  stern  seals  and  a  water  slope 
correction  were  used.   This  last  modification  is  one  which 
has  been  incorporated  into  the  XR-3  program  by  LCDR  G.  T. 
Forbes,  USN  during  his  concurrent  study  and  is  reported  in 
Ref.  5.   The  water  slope  modification  is  used  to  correct 
the  water  level  in  the  plenum  chamber  based  upon  the  fact 
that  increased  pressure  in  the  plenum  will  tend  to  displace 
water  after  it  comes  in  contact  with  the  bow  seal  and  passes 
into  the  chamber.   This  displacement  since  it  does  not  take 
place  instantaneously  is  a  function  of  both  time  and  pressure 
and  hence  results  in  a  "slope"  in  the  water  as  it  passes 
through  the  plenum  chamber.   The  purpose  of  using  both  the' 
original  and  the  new  models  was  to  gain  an  appreciation  of  the 
effects  of  the  modeling  changes  made. 

During  all  the  runs  made  at  Lake  San  Antonio  with  the 
XR-3  the  v;ater  conditions  were  fairly  calm,  but  the  steady 
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state  roll  and  pitch  oscillations  were  always  present  in 
the  processed  data  as  reported  previously.   Figures  19 
through  35  which  are  the  graphs  for  all  the  runs  clearly 
show  this  oscillation  in  both  pitch  and  roll.   Again  it  is 
emphasized  that  these  oscillations  precluded  exact  measure- 
ments of  the  decay  time  constants. 

The  method  used  to  induce  roll  and  pitch  transient 
oscillations  in  the  computer  models  was  to  offset  the  initial 
roll  and  pitch  angle  from  the  known  steady  state  angles 
when  the  program  was  executed.   These  offsets  would  then 
decay  as  the  craft  returned  to  the  steady  state  conditions 
for  that  speed.   The  magnitude  of  these  offsets  were  taken 
from  the  recorded  data.   The  point  in  time  where  this  data 
was  taken  is  marked  on  each  graph  of  the  test  craft's 
recorded  data  by  a  dashed  vertical  line. 

Prom  known  loading  conditions  the  craft's  weight  during 
these  runs  was  605O  lbs.   Its  center  of  gravity  was  located 
10.08  feet  forward  of  the  transJm.   Average  plenum  pressure 
was  not  available  as  this  data  was  not  recorded.   The  plenum 
pressure  was  assumed  to  be  the  nominal  2^1.86  PSF  as  reported 
in  Ref.  3. 

These  inputs  in  the  case  of  the  original  [Ref.  3]  XR-3 
model  caused  the  craft  to  have  a  negative  steady  state  trim 
angle  which  is  readily  apparent  ivhen  looking  at  the  graphs 
of  the  simulation  runs  of  the  original  model  (see  Figures 
22  through  28  for  simulation  runs  at  15-93  knots). 
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Although  it  was  recognized  that  this  negative  trim  angle  was 
not  realistic  the  Initial  conditions  were  not  varied  to  get 
a  positive  trim  angle  because  It  was  desired  to  test  each 
model  with  exactly  the  same  physical  characteristics.   The 
actual  value  of  the  steady  state  trim  was  not  the  quantity 
that  was  to  be  studied.   The  quantity  that  was  of  Interest 
was  the  transient  response  of  the  model  as  It  arrived  at 
that  steady  state  value  after  an  Initial  offset. 

B.   RUN  DATA  ANALYSIS 

The  first  run  v;as  a  calm  water  transient  roll  motion 
study  at  11.79  knots.   In  this  run  the  pitch  offset  was 
considered  negligible.   The  roll  induced  oscillation  was 
measured  to  have  a  negative  magnitude  of  approximately 
.5  degrees.   This  roll  offset  was  entered  into  both  models 
with  the  following  results.   The  decay  time  constant  for 
roll  motion  was  4.6  seconds  for  the  original  XR-3  model  and 
5.6  seconds  for  the  new  model.   The  oscillation  frequency 
was  .752  Hz  for  the  original  XR-3  and  .714  Hz  for  the  new. 
Measured  data  gave  the  actual  boat  a  frequency  of  roll 
oscillation  of  .29  Hz  and  a  decay  time  constant  of  less 
than  1  second. 

The  second  run  was  at  15-93  knots  in  calm  water.   This 
run  had  both  roll  and  pitch  perturbations.   They  were 
minus  .25  degrees  for  pitch  and  a  negative  .27  degrees  for 
roll.   With  these  offsets  the  new  XR-3  model  gave  a  roll 
time  constant  of  4.1  seconds  and  a.  pitch  time  constant  of 
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3.^  seconds.   The  frequencies  of  oscillation  for  the  new 
model  were  .72  Hz  for  roll  and  .4^1  Hz  for  pitch.   These 
same  quantities  were  measured  on  the  old  model  to  be 
3.1  seconds  and  5.0  seconds  for  the  roll  and  pitch  time 
constants  respectively,  and  In  order,  .JS   Hz  and  .50  Hz  for 
the  roll  and  pitch  frequencies.   The  same  measured  quantities 
from  the  recorded  data  gave  something  less  than  one  second 
for  both  the  roll  and  pitch  time  constants  and  .25  Hz  and 
,22  Hz  for  the  oscillation  frequencies  for  roll  and  pitch. 

The  last  calm  water  data  was  taken  at  19.^5  knots.   The 
measured  pitch  offset  was  a  positive  .12  degree  and  a 
positive  .^5  degree  for  the  roll  offset.   Again  these 
values  were  imposed  upon  the  initial  steady  state  trim 
conditions  to  achieve  a  transient  response. 

Prom  the  two  computer  simulations  the  following  results 
were  obtained.   Original  model  roll  frequency  —  .16   Hz, 
roll  decay  time  constant  —  2.7  seconds,  pitch  frequency  — 
.50  Hz,  pitch  time  constant  —  4.6  seconds.   New  model  roll 
frequency  is  .72  Hz  and  decay  time  constant  4.3  seconds. 
Pitch  frequency  is  . 42  Hz  and  decay  time  constant  is 
4.4  seconds.   The  same  data  from  the  recorded  tapes  gave  a 
pitch  frequency  of  .28  Hz  and  with  decay  time  constant  of 
less  than  one  second.   Roll  frequency  was  measured  to  be 
.28  Hz,  and  the  roll  decay  time  constant  was  estimated  at 
less  than  one  second. 
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From  the  above  data  taken  on  the  three  runs  and  the 
computer  simulations  thereof  It  is  evident  that  the  computer 
model  neither  in  the  original  form  nor  in  its  new  form 
accurately  predicts  roll  or  pitch  transient  motion  under 
the  conditions  studied.   Decay  time  constants  for  pitch 
and  roll  were  generally  four  times  too  long  in  both  simula- 
tion models.   Pitch  and  roll  frequencies  of  the  computer 
models  were  generally  too  high  by  a  factor  of  two.   Table 
III  gives  this  transient  behavior  of  the  XR-3  and  both 
simulation  models  in  tabular  form. 

Figure  19  is  the  actual  roll  response  data  for  Run  1 
collected  at  11.79  knots.   Figures  20  and  21  are  the 
computer  simulation  runs  for  the  original  and  new  model 
for  this  run  depicting  roll  response.   Recorded  pitch  and 
roll  motion  for  Run  2  (15.93  knots)  is  given  in  Figures  22, 
23  and  2^  and  the  original  model  computer  simulation  of 
roll  and  pitch  Is  in  Figures  25  and  26  respectively.   The 
response  of  the  new  model  for  Run  2  is  recorded  in  graphs 
27  and  28.   Recorded  data  for  Run  3  Is  given  in  Figures  29 
through  31  with  the  computer  model's  simulations  of  pitch 
and  roll  given  In  Figures  32  through  35. 

C.   SIMULATION  STUDIES  OF  THE  TWO  MODELS 

In  an  effort  to  gain  better  knowledge  of  the  total 
effect  of  the  new  modeling,  both  the  old  and  the  new  models 
were  given  a  positive  three  degree  pitch  offset  from  steady 
state  calm  water  conditions.   Both  models  were  set  at  the 
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TABLE    III.       CALM  WATER  FORCED  OSCILLATIONS   RESPONSE 


Run   1,    Speed   11.79   knots 


Variable 

Roll  Frequency 

Roll  decay  tine  constant 

Pitch  Frequency 

Pitch  decay  time  constant 


Measured 
XR3 

.29  Hz 
<1  sec 


Simulated 
Original  MDdel   New  Model 


.752  Hz 
4.6  sec 


.714  Hz 
5.6  sec 


Run   2,   Speed   15.93  knots 


Variable 

Roll  Frequency 

Roll  decay  time  constant 

Pitch  Frequency 

Pitch  decay  time  constant 


rfeasured 
XR3 

.25  Hz 

<1  sec 

.22  Hz 
<1  sec 


Simulated 
Original  Model      New  Model 


.76  Hz 
3.1  sec 
.50  Hz 
5.0  sec 


.72  Hz 

4.1  sec 
.44  Hz 
3.4 


Run    3,   Speed  19.45  knots 


Variable 

Roll  Frequency 

Roll  decay  tiire  constant 

Pitch  Frequency 

Pitch  decay  time  constant 


ffeasured 
XR3 

.28  Hz 
<1  sec 

.28  Hz 
<1  sec 


Simulated 
Original  Model   New  Model 


.76  Hz 
2.7  sec 
.50  Hz 
4.6  sec 


.72  Hz 

4.3  sec 
.42  Hz 

4.4  sec 
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FIGURE   19 

Roll  Transient  verses  Time^  Velocity  11.79  knots 
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FIGURE  20 

Computer  Simulation  Roll  Angle  verses  Time,  11.79  knots 
Original  XR-3 
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FIGURE    21 

Computer  Simulation  Roll  Angle  verses  Time.  11.79  knots 
New  XR-3 
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FIGURE   22 

Pitch  verses  Time,  Velocity  15.93  knots 
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FIGURE  23'  . 

Pitch  Transient  verses  Time,  Velocity  15.93  knots 
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FIGURE    2H 

Roll  Transient  verses  Time^  Velocity  15.93  knots 
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FIGURE  25 

Computer  Simulation  Pitch  Angle  verses  Time^  15.93  knots 
Original  XR-3 
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FIGURE   26  ' 

Computer  Simulation  Roll  Angle  verses  Time,  15.93  knots 
Original  XR-3 
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FIGURE  2  7 

Computer  Simulation  Pitch  Angle  verses  Time,  15.93  knots 
New  XR-3 
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FIGURE  2  8 

Computer  Simulation  Roll  Angle  verses  Time, 15. 93  knots 
New  XR-3 
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FIGURE   2  9 

Roll  Transient  verses  Time^  Velocity  19.45  knots 
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FIGURE    30 

Pitch  verses  Time',  Velocity  19.^5  knots 
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FIGURE  31 

Pitch  Transient  verses  Time.  Velocity  19.^5  knots 
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FIGURE  32 

Computer  Simulation  Roll  Angle  verses  Time,  19.^5  knots 
Original  XR-3 
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FIGURE   33 

Computer  Simulation  Pitch  Angle  verses  Time^  19.^5  knots 
Original  XR-3 


92 


Roll 
Angle 
(degrees' 
(x.l) 


-?  ^ 


FIGURE  3^ 

Computer  Simulation  Roll  Angle  verses  Time,  19.^5  knots 
New  XR-3 
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FIGURE  35 

Computer  Simulation  Pitch  Angle  verses  Time,  19.45  knots 
New  XR-3 
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same  initial  conditions  that  they  were  during  the  19.45  knot 
run  except  for  the  offset.   In  the  case  of  the  new  model 
this  raised  the  bow  approximately  six  Inches  and  since  the 
steady  state  condition  gave  the  bow  a  depth  of  about  three 
Inches  with  the  positive  steady  state  pitch  angle  this  caused 
air  gaps  to  appear  beneath  the  bow  seal.   In  the  case  of  the 
old  model  where  the  steady  state  pitch  angle  was  negative 
no.  air  gaps  were  created.   Figures  36  through  38  present  the 
original  model's  response  and  Figures  39  through  4l  present 
the  new  model's  response. 

As  can  be  seen  in  Figure  40,  the  air  gap  on  the  new  model 
caused  wide  fluctuations  in  the  plenum  pressure  during  the 
first  half-second.   This  loss  of  pressure  permitted  the  boat 
to  rapidly  sink  as  shown  in  Figure  hi   until  the  plenum 
pressure  built  up.   When  the  air  gaps  closed  and  the  boat 
continued  to  settle  the  plenum  pressure  rapidly  rose  as  the 
plenum  volume  decreased.   The  boat  exponentially  returned 
to  steady  state  draft  levels  as  the  pitch  motion  shown  in 
Figure  39  died  down.   The  pitch  motion  of  the  new  model 
shows  the  effect  of  the  new  seals.   The  downward  movement  . 
of  the  bow  as  seen  in  Figure  39  is  severely  limited  as  the 
bov7  seal  wetted  length  becomes  greater  and  the  effect  of 
the  rigid  body  analysis  comes  into  play.   The  first  negative 
overshoot  is  limited  to  34  percent  of  the  original  effect. 
The  initial  decay  time  constant  Is  2.5  seconds  and  the 
pitch  frequency  is  .63  Hz  which  because  of  the  nonlinearlty 
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FIGURE  36 

Computer  Simulation  3  Degree  Pitch  Step  Response^  19.^5  knots 
Original  XR-3   Pitch  verses  Time 
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FIGURE    37 

Computer  Simulation  3  Degree  Pitch  Step  Response.  19.45  knots 
Original  XR-3  Plenum  Pressure  verses.  Time 
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FIGURE  38 

Computer  Simulation  3  Degree  Pitch  Step  Response^  19.^5  knots 
Original  XR-3  Draft  at  C.G.  verses  Time 
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FIGURE  39 

Computer  Simulation  3  Degree  Pitch  Step  Response,  19.^5  knots 
New  XR-3  Pitch  verses  Time 
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FIGURE  ^0 

Computer  Simulation  3  Degree  Pitch  Step  Response,  19.^5  knots 
New  XR-3  Plenum  Pressure  verses  Time. 
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FIGURE  41 

Computer  Simulation  3  Degree  Pitch  Step  Response^  19.^5  knots 
New  XR-3  Draft  at  C.6.  verses  Time 
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of  the  bow  and  stern  seal  models  changes  to  a  frequency  of 
.k2   Hz  at  the  end  of  ten  seconds.   The  time  constant  has 
lengthened  to  4.3  seconds  In  the  same  period. 

The  behavior  of  the  old  model  is  different.   Because 
of  the  high  capacity  of  the  fan  maps  the  damping  effect  of 
plenum  pressure  is  minimal.   The  draft  movement  shown  in 
Figure  38  is  therefore  highly  underdamped  whereas  the  draft 
motion  for  the  new  model  did  not  display  any  oscillatory 
tendencies  (see  Figure  41)  and  displayed  exponential  decay 
type  response  when  returning  to  steady  state  conditions. 
In  pitch  motion  shown  in  Figure  36  the  original  model 
displays  a  first  overshoot  magnitude  of  78  percent  of 
original  offset.   The  frequency  of  oscillation  is  .5  Hz 
and  the  decay  time  constant  is  4.78  seconds.   Both  the  time 
constant  and  pitch  oscillation  frequency  remain  constant 
throughout  the  IC  second  simulation. 

To  further  evaluate  the  differences  between  the  original 
and  new  models  a  second  simulation  run  was  made  at  19-45 
knots.   This  run  was  in  ahead  sea  state  2  irregular  seas. 
Figures  42  and  43  show  that  the  plenum  pressure  fluctuations 
in  the  new  model  were  6OO  percent  greater  than  in  the 
original  model.   VJhereas  the  plenum  pressure  varied  from 
5  to  about  45  PSF  in  the  new  model  the  pressure  only  varied 
from  about  21  to  27  PSF  in  the  original.   Draft  fluctuations 
were  about  300  percent  greater  in  the  new  model  than  in  the 
original  (see  Figures  44  and  45).   Pitch  angle  fluctuations 
were  comparable  in  both  models  (Figures  46  and  47). 
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FIGURE    il2 

Computer  Simulation  Plenum  Pressure  verses  Time,  19. '45  knots 
New  XR-3;  Sea  State  2 
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FIGURE    i<3 

Computer  Simulation  Plenum  Pressure  verses  Time^  19.^15  knots 
Original  XR-3.  Sea  State  2 
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FIGURE  ^i< 

Computer  Simulation  Draft  at  C.G.  verses  Time.  19.45  knots 
New  XR--3;  Sea  State  2 


105 


Draft 
(inches ) 


I  to  tz 

Time     (seconds) 


ft 


FIGURE    i|5 

Computer  Simulation  Draft  at  C.G.  verses  Time.  19.^5 
Original  XR-3,  Sea  State  2 
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FIGURE    hS  ■ 

Computer  Simulation  Pitch  Angle  verses  Time,  19.^5  knots 
New  XR-3;  Sea  State  2 
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FIGURE  ^1 

Computer  Simulation  Pitch  Angle  verses  Time^  19.^5  knots 
Original  XR--3,  Sea  State  2 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

Making  accurate  measurements  from  the  collected  data 
was  extremely  difficult  because  of  the  high  relative  level 
of  steady  state  oscillation  present  In  both  the  pitch  and 
roll  motion.   Signal  noise  coupled  with  the  fact  that  the 
decay  time  constants  and  the  periods  of  the  oscillations 
were  of  the  same  general  magnitude  made  accurate  measurement 
of  the  decay  time  constants  practically  Impossible. 

Although  accurate  measurement  of  the  time  constants  was 
not  possible  general  observations  could  be  made.   It  was 
clearly  evident  that  the  pitch  and  roll  transients  Induced 
In  the  XR-3  test  craft  decayed  very  rapidly  and  not  with 
the  highly  underdamped  characteristics  shown  by  the  computer 
simulation  models.   The  measured  tape  data  confirmed  the 
observations  made  aboard  the  XR-3  while  conducting  these 
runs  that  the  craft  come  to  steady  state  conditions  with 
only  one  or  two  overshoots  after  the  forcing  action  had 
ceased.   Roll  motion  transients  were  normally  completed  in 
less  than  two  cycles.   Both  the  longitudinal  and  athwartshlps 
stability  were  highly  speed  dependent  gaining  stability  as 
speed  increased. 

These  observations  were  not  based  solely  on  "feel"  but 
were  made  by  sighting  across  the  deck  of  the  XR-3  and 
observing  the  motion  of  the  deck  edge  relative  to  the  shore 
line. 
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Simulation  runs  have  shown  that  for  large  perturbations 
the  transient  response  of  the  new  model  Is  quite  different 
from  the  original.   Since  It  was  not  possible  to  excite  such 
large  perturbations  in  the  actual  test  craft  It  remains  to 
be  verified  that  the  new  bow  and  stern  seal  models  more 
accurately  predict  the  forces  generated  by  the  actual  seals 
than  does  the  original  seal  model. 

The  new  seal  fan  maps  are  more  accurate  than  the  original 
since  they  are  based  on  actual  verified  manufacturer's  data. 
Leakage  areas  are  now  proportional  to  those  of  the  actual 
craft. 

The  fact  that  the  present  computer  models  of  the  XR-3 
both  old  and  new  do  not  appear  to  accurately  represent  the 
actual  transient  behavior  of  the  XR-3  in  pitch  and  roll 
opens  many  areas  which  should  be  investigated  for  possible 
errors.   As  both  the  pitch  and  roll  transients  appear  to  be 
in  error  the  sidewall  subroutine  becomes  a  prime  area  of 
investigation  as  it  contributes  major  components  to  both 
the  roll  and  pitch  force  moments. 

Figure  ^8  shows  one  area  which  should  be  explored  and 
that  is  the  area  of  hydrodynamic  lift  of  the  sidewalls. 
Figure  ^S   is  a  simplified  drawing  of  the  underside  of 
the  XR-3.    The  flat  sections  on  the  bottom  of  each 
sidewall  are  cross-hatched  for  clarity.   At  high 
velocity  (20-25   knots)  these  areas  become  effective 
planing  surfaces.   In  this  program  the  lift  of  the 
sidewall  sections  due  to  planing  does  not  appear  to  take 
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Into  account  the  underwater  body  shape.   This  shape  changes 
radically  from  bow  to  stern.   Whether  or  not  this  shape  has 
significant  effect  on  roll  and  pitch  transient  behavior  is 
not  known  but  should  be  investigated.   Another  aspect  of 
the  sidev/all  modeling  which  should  be  investigated  is  the 
Inclusion  of  forces  due  to  the  displacement  of  water  per 
unit  of  time.   A  term  such  as  this  would  not  only  have  terms 
based  on  draft  and  pitch  angle  but  would  also  include  terms 
which  took  into  account  the  rate  of  change  of  displacement 
and  the  rate  of  change  of  the  pitch  angle.   The  addition  of 
first  derivative  terms  such  as  these  would  have  a  definite 
effect  on  the  transient  behavior  of  the  sidewalls. 

A  third  area  of  study  concerning  the  sidewalls  is  the 
effect  of  sidewall  deadrise  on  turning  performance.   Modeling 
of  the  deadrise  forces  would  be  incorporated  in  the  cross- 
flow  force  modeling.   This  deadrise  could  have  significant 
effect  on  the  outboard  sidewall  during  a  turn  determining 
whether  it  has  a  tendency  to  rise  and  ride  over  the  water 
as  the  craft  skids  or  dig  in  and  cause  the  craft  to  undergo 
violent  roll  motion  as  reported  in  Ref.  6. 

Another  area  of  study  which  needs  investigation  is  the 
effect  of  plenum  pressure  on  the  transient  response  (of  the 
computer  model).   If  it  can  be  shown  that  the  computer 
model  as  presently  configured  can  be  made  to  duplicate 
actual  measured  craft  transient  motion  by  Just  adjusting 
the  plenum  pressure,  further  modifications  of  the  program 
may  be  unnecessary.   Further  transient  motion  studies  are 
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needed  where  both  the  plenum  pressure  and  the  seal  pressures 
are  Included  in  the  recorded  data.   Because  of  the  stochastic 
nature  of  the  actual  measured  data  many  duplicate  runs  will 
be  needed  to  arrive  at  accurate  averages  of  the  test  craft's 
oscillation  frequencies  and  the  decay  time  constants.   These 
averages  should  then  be  used  in  comparing  the  actual  craft 
with  the  computer  model. 

As  can  be  seen  in  the  recorded  data  presented  herein  a 
fruitful  area  of  study  would  be  in  reducing  the  signal  noise 
on  the  recorded  data.   This  could  be  in  the  recording  phase 
or  processing  phase  of  the  data.   Thorough  tests  under 
controlled  conditions  need  to  be  made  of  the  actual  transfer 
characteristics  of  the  different  sensors  installed  in  the 
test  craft. 

Earlier  in  this  report  it  was  mentioned  that  there 
existed  steady  state  roll  and  pitch  oscillations  while  the 
boat  was  traveling  straight  on  calm  water.   These  oscillations 
could  exist  and  could  be  the  result  of  the  craft  being  in 
some  limit  cycle  type  of  motion  because  of  the  inherent 
nonlinearity  of  vehicles  operating  in  fluid  mediums.   On 
the  other  hand,  it  could  be  the  result  of  instrumentation 
for  example  if  the  vertical  seeking  gyro  from  which  all  roll 
and  pitch  motions  are  derived  was  precessing  about  the  true 
vertical  then  it  could  impose  a  sinusoidal  oscillation 
signal  on  both  the  pitch  and  roll  sensors  that  would  show 
up  in  the  recorded  data. 


ii: 


In  conclusion  It  must  be  stated  that  there  are  still 
many  areas  in  this  subject  of  sidewall  modeling  of  the 
XR-3  which  should  and  need  to  be  explored  before  the  XR-3 
model  as  originally  stated  in  Ref.  3  and  modified  herein 
can  be  accepted  as  a  reasonably  accurate  representation  of 
the  actual  XR-3  test  craft.   Further  investigation  into 
the  data  collection  and  processing  is  also  needed  to  refine 
the  processed  data  to  the  point  where  it  is  accurate 
enough  for  dynamic  studies. 
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APPENDIX  A 
LISTING  OF  BOWSEAL  SUBROUTINE  'bOWSl' 


The  entering  arguments  needed  for  the  new  bow  seal 
routine  are  the  same  as  presented  In  Ref.  3  for  the  original 
bow  seal  except  that  Instead  of  entering  the  angle  between 
the  leading  edge  and  craft  vertical  in  columns  46-55,  the 
length  of  the  rear  cables  Is  entered  (In  feet).   An  additional 
card  Is  needed  to  enter  the  length  of  the  middle  support 
cables.   The  new  data  card  should  be  entered  as  shown  below. 


CARD 


COLUMNS 

FORMAT 

1-  5 

I 

6-15 

P 

16-25 

P 

26-35 

P 

36-^5 

P 

H6-55 

56-65 

66-75 

1-  5 

6-15 

ENTRY 

00601  -  control  tag 

X  coordinate  of  seal 
hinge  (ft)  forward 
of  transom 

Seal  leakage  orifice 
coefficient 

Pressure  differential 
between  bow  seal  and 
plenum  chamber  (PSF) 

Z  coordinate  of  seal 
hinge  (ft)  above  keel 

Length  of  rear  support 

cable  (ft) 

(1.875  is  full  down) 

Length  of  the  seal 
leading  edge  (ft) 

Base  leakage  area  (ft  ) 
Nominal  .08  ft^) 

00602  -  control  tag 

Length  of  middle  support 

cable  (ft) 

(1.1  R7f^  1s  full  down) 
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APPENDIX  B 
LISTING  OF  STERN  SEAL  SUBROUTINE  STNSL 

The  entries  for  this  routine  are  the  same  as  the  original 
stern  seal  except  that  in  columns  46-55  the  rear  support 
cable  length  is  entered  instead  of  the  seal  angle. 

Block  5  entries  are  as  follows: 

ENTRY 

00500  -  control  tag 

X  coordinate  of  seal  hinge, 
(ft)  forward  of  transom 

Z  coordinate  of  seal  hinge 
above  keel  (ft) 

Base  leakage  area  (ft  ) 
3.79  nominal 

Seal  leakage  orifice 
coefficient 

Length  of  rear  support 
cable  (ft)  1.875  is  full 
down 

56-65         F         Pressure  differential 

between  stern  seal  and 
plenum  chamber  (PSF) 

66-75         F         Length  of  the  leading  edge 

of  the  seal  (ft) 


CARD 

COLUMNS 

FORMAT 

1 

1-  5 

I 

6-15 

F 

16-25 

F 

26-35 

F 

36-45 

F 

46-55 

F 
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APPENDIX  B 
LISTING  OP  STERN  SEAL  SUBROUTINE  STNSL 
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APPENDIX  C 
LISTING  OF  FAN  SUBROUTINE  PAN* 

The  entries  for  the  fan  subroutine  remain  the  same. 
The  nominal  fan  RPM  is  8000  RPM  vice  2615  RPM.  Printed 
below  are  copies  of  the  appropriate  entries  for  the  bow 
seal,  stern  seal  and  main  plenum  fans. 


01901   2. 
-10.406 
23.414 

0 

8000.0   19. 
-5.203     0.0 
26.015     28.617 

0 

1.0 
5.203 
31.218 

10.406 
33.820 

15.609 
36.421 

18.211 
41.624 

20.812 
46.82  7 

52.030 

22.7 

12.20 

57.233 
22.067 
10.00 

62.436 

20.67 

7.80 

19.60 
6.40 

18.367 
5.23 

17.133 
4.33 

16.467 
2.60 

15.10 
.983 

-.020 

-1.06 

-2.03 

01902   2. 
-10.406 
23.414 

.0 

8000.0   19, 
-5.203     0.0 
26.015     28.617 

,0 

1.0 
5.203 
31.218 

10.406 
33.820 

15.609 
36.421 

18.211 
41.624 

20.812 
46.827 

52.030 

22.7 

12.20 

57.233 
22.067 
10.00 

62.43b 

20.67 

7.80 

19.60 
6.40 

18.367 
5.20 

17.133 
4.33 

16.467 
2.60 

15.10 
.983 

-.020 

-1.06 

-2.03 

01903   1. 
-10.406 
23.414 

.0 

8000 
-5.203 
26.015 

.0    19. 
0.0 
28.617 

,0 

1.0 
5.203 
31.218 

10.406 
33.320 

15.609 
36.421 

18.211 
41.624 

20.812 
46.827 

52.030 
25.70 
16.40 
0.00 

57.233 
24.33 
1Z.60 
-1.06 

62.436 
23.167 
10.367 
-2.03 

21.83 
8.033 

20.61 
6.331 

"^  19.33 
5.073 

13.85 
3.035 

17.73 
1.60 
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APPENDIX   C 
LISTING   OF   PAN   SUBROUTINE   FAN 
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APPENDIX  D 
LISTING  OF  INITIAL  CONDITION  SUBROUTINE  INCOn' 


The  INCON  subroutine  was  modified  to  permit  it  to 
calculate  the  wave  components  of  the  various  sea  states  if 
desired  instead  of  having  to  enter  each  component  individually. 

Block  11  options  1  and  2  remain  identical.   New  options 
3  and  H   have  been  added  and  are  entered  as  follows: 
BLOCK  11  -  OPTION  TAG  3 


CARD 


COLUMS 

FORMAT 

1-   5 

I 

6-15 

F 

16-25 

P 

26-35 

F 

36-^5 

F 

H6-55 

F 

ENTRY 

01103  control  and  option  tag 

Number  of  wave  components 
maximum  of  10 

Initial  wave  heading  (deg) 
(l80  deg  =  ahead  waves) 

Average  height  of  1/3 
highest  waves  (ft) 

Shortest  significant  period 
(sec) 

Longest  significant  period 
(sec) 


BLOCK    11    -    OPTION 

TAG    4 

CARD              COLUMNS 

FORMAT 

1                   1-5 

I 

6-15 

F 

16-25 

F 

26-35 

F 

ENTRY 

01104  control  and  option  tag 

Number  of  wave  components 
maximum  of  10 

Initial  wave  heading  (deg) 
(l80  deg  =  ahead  vjaves) 

Average  height  of  1/3 
highest  v;aves  (ft) 
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CARD 

COLUMNS 

FORMAT 

ENTRY 

1 

36-45 

P 

Lowest  significant  frequency 
(Hz) 

116-55 

F 

Highest  significant  frequency 
(Hz) 
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APPENDIX  D 
LISTING   OF   INITIAL   CONDITION   SUBROUTINE   INCON 
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